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EXPERIMENTAL AND ANALYTICAL STUDIES OF 
A TRUE AIRSPEED SENSOR 

By 

G.L. Goglia^ cuid J.Y. Shen^ 
ABSTRACT 


The objective of this study was to analyze and design a true airspeed 
sensor to replace the conventional pitot-static pressure transducer. 

This indicator should have a flow phenomenon which is the vortex preces- 
sion or the "vortex whistle," no moving parts, and also be independent 
of temperature, density, altitude or h\amidity changes- This indicator 
was to be designed mainly for small commercial aircraft with airspeeds up 
to 321.9 km/hr (200 mph) . 

In an attempt to model the complicated fluid mechanics of the vortex 
precession, three-dimensional, inviscid, unsteady, incompressible fluid 
flow was studied by using the hydrodynamical linearized stability theory. 
The temporal stcibility approach was used to derive the relationship 
between the true airspeed cind frequency response- The results show that 
the frequency response is linearly proportional to the airspeed. 

The designed sensor basically consisted of a vortex tube, a swirler, 
and a pickup system. When air passed through the swirler, a precessional 
flow was generated at the region before and after the sudden enlargement 
area. An audible vortex whistle was also generated. The signal was 
picked up by the microphone, auid the frequency response shown in a 
frequency counter. Measurements for both the closed-conduit test and wind 
tunnel test were recorded. 

A computer progr 2 un was developed to obtain the numerical solution. 

The parameters describing the sensor were introduced into the calculation. 


Eminent Professor eind Chairperson, Department of Mechanical Engineering 
and Mechcuiics, Old Dominion University, Norfolk, Virginia 23508 

^ Graduate Research Assistant, Department of Mechanical Engineering and 
Mechanics, Old Dominion University, Norfolk, Virginia 23508 


Computational results for various combinations of vortex tubes and swirlers 
were obtained. 

For a specific flow rate or airspeed, the larger the exit swirler 
angle, the greater was the frequency response. For a smaller cross- 
sectional area at a precessional flow region, the frequency response was 
higher. It was observed that as the airspeed was increased the Strouhal 
number remained constant. The Strouhal number was found to be only 
dependent on the exit angle of swirler. 

In some cases, the experimental results were found to be in reasonable 
agreement with the theoretical predictions. The viscous effect on the 
performance of the sensor was apprecicdsle , especially for some sensors. 
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LIST OP SYMBOLS 


,C 


Cl , C2 


Di 

D2 

D3 

f 

F 

g 

h 

k 

K 

n 

P 

P' 


I 


I 


q 

Q 

r 

R 

Re 

St 

t 


constant, defined in equation (52) 
arbitrary constcuits 

diameter of swirler generation region, defined in figvire 1 , cm 

diameter- of swirler, defined in figure 1 , cm 

diameter of processional flow region, defined in figure 1 , cm 

function of radius r, amplitude of preliminary complex solution 
in r-direction 

frequency response, KHz 

function of radius r, amplitude of preliminary complex solution 
in 0 -direction 

function of radius r, amplitude of preliminary complex solution 
in z-direction 

wavenumber, a real number 

parameter, defined in equation (40) 

integer 

actual fluid pressure, kg/cm^ 
pressure fluctuation, kg/cm^ 
preliminary complex solution in pressure 
mean pressure, kg/cm^ 

function of radius r, amplitude of prelimineury complex solution 
in pressure 

flow rate, m^/min 

radius or radial coordinate, cm 

radius of precessional flow region, defined in figvure 1 , cm 
Reynolds nximber based on mean axial velocity and Da 
Strouhal number based on mean axial velocity and D 3 
time, sec 
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actual fluid radial velocity, km/hr 
velocity fluctuation in r-direction, km/hr 
preliminary complex solution in r-direction, km/hr 
true airspeed, km/hr 

actual fluid tangential velocity, km/hr 
velocity fluctuation in 6-direction, km/hr 
preliminary complex solution in 0-direction, km/hr 
mean velocity in 0-direction, km/hr 
actual fluid axial velocity, km/hr 
velocity fluctuation in z-direction, km/hr 
preliminary complex solution in z-direction, km/hr 
mean velocity in z-direction, km/hr 
constant mecin axial velocity, km/hr 
axial coordinate 

exit swirler angle, defined in figure 59 

6-coordinate 

fluid density, g/cm^ 

angular velocity, a constant complex number 
constant upstream fluid angular velocity 
complex conjugate of preliminary complex solution 



INTRODUCTION 


At the present time there are two types of airspeed instruments which 
might be used on aircraft: the differential-pressure type and the true 

airspeed meter. However, the pitot-static instrument, which is of the 
differential-pressure type, is exclusively used. This indicator is 
calibrated in terms of airspeed at a standard air density. In order to 
obtain the actual airspeed at other densities, a correction must be 
made. The performance of the pitot-static tube is greatly affected by 
installation location. One must therefore find a location for the pitot- 
static openings that will be free from structual interference effects. 

The true airspeed sensor is the conventional type of meter with 
rotating surfaces, such as propellers, which gives readings independent of 
air density. The sensor is usually used in making measurements of air- 
speeds in the lower ranges. One kind of true airspeed sensor used by 
the United States Navy on airships is known as the commutator-condensor 
type. 

The idea of designing a true airspeed sensor originated from the 
discovery of vortex whistle and the flow phenomenon— precession , which 
is different from that of vortex shedding. The origin of the sound or 
whistle was studied by a few investigators (refs. 1-3) from different 
points of view. The vortex shedding problems associated with the 
aerolian tones and edgetones were also studied aind observed experimentally 
by many researchers (refs. 4-12) . Very few people worked on the vortex 
whistle and precessional flow problems. Bernard Vonnegut (ref. 13) in 
1954 was the first to discover and investigate the vortex whistle. 

His laboratory was conducting an experiment on a vortex creating housing 
for aircraft thermometers. During the experiment he observed that a 
sound was generated when the rotating air escaped from the open end of 
the tube. He also found that the frequency of this sound increased 
with increasing rates of air flow. In addition, the frequency that was 
produced decreased as the length of the tube in which the vortex rotated 
was increased. Vonnegut suggested that the vortex instability leaving 
the tube caused the whistle and developed an empirical formula describing 
the performcince of his whistle. 



In 1955, Irving Michelson (ref. 14) published the first analytical 
work on the theory of a vortex whistle. He considered the flow throughout 
the whistle to be two-dimensional, unsteady, inviscid, and isentropic. 

He introduced the parameters s , u ' , and v ' , where s was a non- 
dimensional density perturbation. By this method he was cible to arrive at 
linearized simultaneous equations. A secular equation was then derived, 
noting one root of particular interest. From the solution and the secular 
equation, he noted the occurrence of a frequency which was proportional 

to the flow speed U. When he introduced the isentropic flow relationship, 
he was able to express the frequency in terms of pressure drop and 
reservoir sound speed Cj . 

Michelson 's theory compares very favorably with Vonnegut's empirical 
formula, particulairly when constants 1.09 and 1.19 are used in Michelson 's 
theory for monoatomic and diatomic gases, respectively. 

In 1957, J.P. Nicklas (ref. 15) reported his work on the investigation 
of a vortex tube, acoustic, tirue airspeed sensor conducted at the Cornell 
Aeronautical Ladsoratory. He investigated the feasibility of measuring 
true aircraft airspeed by measuring the frequency of the sound produced 
in a vortex tube mounted on an airplane. Nicklas, however, concentrated 
his efforts on the single tangential nozzle vortex tube. His data 
revealed that the fundamental sound frequency of a vortex tube could be 
considered a linear function of true airspeed in the subsonic speed 
range. He indicated that the altitude and temperature sensitivities 
of the vortex tube could be reduced by proper design. In his conclusions, 
he also mentioned that, no significant improvement in signal quality was 
obtained by modifying the txibe shape. Nicklas also studied the effect 
the angle of attack had on frequency response. 

In 1960, M. Suzuki (ref. 16) conducted studies to determine a method 
of eliminating the whistle occurring in the vortex tube. In his analysis 
he assumed both a free and forced vortex region of velocity distribution. 
Suzuki, using the boundary conditions at the wall and at the interface 
between free and forced vortex, derived a linearly proportional relation- 
ship between the peculiar frequency and the angular velocity of the 
forced vortex. In his derivation, the density and velocity components 
were separated into mean and fluctuation terms. Suzuki introduced a 
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number of assumptions and restrictions in order to obtain the Bessel's 
equation and its solution. Although Suzuki did not present either 
numerical or quantitative results, he did report and discuss his experi- 
mental data. In addition to the linear relationship, Suzuki found that 
no soxind was produced at small flow rates and that, when the value of 

L /D was less than unity, no distinct frequency could be observed; L 
c c '' 

was the length of cold tube in his model, and was the "diameter 

of the outlet. 

In 1963, Robert C. Chanaud (ref. 17) converted Vonnegut's data 
into Reynolds and Strovihal numbers and found that the air euid water data 
were almost coincident, suggesting that dynamic similarity might occur. 

The perturbation of a two-dimensional inviscid vortex flow was investigated. 
Chanaud derived a linear relationship between pertvirbation frequency and 
fluid angular velocity for neutrally stable oscillations of an inviscid 
flow. His results support the investigations and conclusions reported 
by Vonnegut and Michelson. He confirmed that the precessional frequency 
is the same as the sound frequency and that the fluid angular velocity 
is simply related to the precession frequency of the unstciblc motion. 

In his conclusion, he mentioned that high speed was not necessary to 
generate the whistle, as velocities of 1.5 m (5 ft) per sec were found 
sufficient. He, as others did , explained that the instability which 
occurred was due to the sudden area change at the tube exit. Chanaud 's 
results show that the amplitude of oscillation within the tube depends 
on how the area changes; a gradual area increase permits larger 
amplitude flow oscillations whereas an abrupt area change reduces the 
magnitude of the flow oscillation within the tube. He mentioned 
that this may be the reason Vonnegut did not detect the sound with a 
flared tube, and also stated that no quantitative information on the 
nature of the instability had been obtained. 

Powell (ref. 18) in 1964 pxoblished a paper discussing the origin 
of the sound. He showed and explained in detail, from a physical point 
of view, how aerodynamic sound in an unsteady fluid flow was generated 
as a result of the movement of vortices, or of vorticity. 

In 1965, Chanaud (ref. 19) pxoblished a paper describing the experi- 
mental study in certain swirling flows . One of the swirling flows was 



studied by Talbot (ref. 20). The experimental results show that the 
periodic motion in both a vortex whistle and a cyclone separator can be 
described in terms of a hydrodynamic oscillator where the frequency is 
closely related to the cingular velocity of the flow. Chanaud also 
mentioned that the two important parameters, the Reynolds number and 
the Strouhal number, are both of such magnitude that it appears no 
important simplifications can be made in the equations of motion to solve 
the problem analytically. The energy of the oscillator is derived from 
the hydrodynamic instability of the fluid within a reversed-flow region 
on the swirl axis. No qucintitative information is available on the 
condition of a steady reversed-flow region. Chanaud, however, mentioned 
that the experimental results suggest that the two-dimensional per- 
turbation analysis may prove of some value in describing the amplifier 
part of the oscillator. 

Rodely et al. (ref. 21) found that the oscillative motion began 
only beyond certain Reynolds nxambers. They also observed that the oscil- 
lative motion was accompanied by the reversal of flow near the t\ibe axis. 
Gove and Ranz (ref. 22) explained this reversal of flow in detail. The 
reversal of flow was caused by the sudden area enlargement at the tube 
exit. In the better swirler designs, the Rossby mamber could be held 
constant for various Reynolds numbers. This indicated that the frequency 
was linearly related to the flow rate. However, below some Reynolds 
number, due to viscous effects, there were deviations from the constant 
value . 

Chanaud again in 1970 (ref. 23) suggested that the vibrating system 
in the aerodynamic whistle is the air itself . This is in contrast to 
nonaerodynamic devices such as a drum or loud speaker where sound is 
generated when a mechanical system vibrates and disturbs the air. Chanaud 
showed that, due to the instability of the system, a small disturbance 
in the stream flowing through the aerodynamic whistle was amplified, 
and that kinetic energy was converted to oscillatory energy. Part of 
the energy of the amplified disturbance is fed back upstream where the 
flow is most unstcible, and, if the right frequency and amplitude exists, 
it interacts with the original disturbance to maintain the process. 

g fsw cycles the feedback controls the input completely. A whistle 
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is produced when the flow speed is high enough and the frequency is in 
the audible range. 

As mentioned in references 21 and 24 to 31, there is one common 
feature that introduces the concept of "no moving parts" in fluidic 
devices. In contrast to this concept is the device with moving parts as 
in reference 32. Fluidic devices have been widely researched in the 
past 19 years. Simplicity, reliability and easy maintenance make fluidic 
devices attractive. Reference 33 states, "Although present theory gives 
results sufficiently accurate for engineering design, it is not possible 
to justify all the assumptions used. Thus in a scientific sense the 
theory is not always satisfying, but in an engineering design sense the 
theory does seem to be satisfactory." 

In this investigation fluidic models were designed and then tested 
in both water and air. Flow visualization tests in a water model were 
undertaken in order to actually see the flow phenomenon of precession. 
Smaller models were svibsequently made for testing with compressed air 
and in a wind tunnel. A complete analytical analysis was provided in 
this study. The physical models were simulated and used in computer 
calculation. The numerical solutions involved tnae airspeeds up to 
321.89 km (200 mi) per hr. Six different combinations of vortex tubes 
and swirlers were used both in computer calculations as well as in 
experimental tests. 

The objective of this study was twofold. The first objective was 
to analyze and design a true airspeed indicator to replace the conventional 
pitot-static pressure transducer for small commercial aircraft. The 
desired features of this indicator should include the flow phenomenon- 
precession, vortex whistle, and no moving parts. In addition this 
indicator should not be affected by temperature, density, altitude, and 
humidity changes. The second objective was to obtain a numerical 
solution and predict the frequency response which is generated by the 
vortex whistle at a certain airspeed. Prior to this study, theoretical 
results were never presented quantitatively to enable a comparison with 
experimental data. Thus, this study will present a general solution to 
the problem and also provide specific analytical results for comparison 
purposes. A correction factor for viscous effects will later be 



introduced to enable a correlation between theoretical results and 
experimental data. 
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THEORY 


A three-dimens lonalr Invlscid, unsteady, incompressible flow was 
studied by using the hydrodynamical linearized stability theory. The 
hydrodynamical linearized stcdaility theory, as mentioned in references 
34 to 40, has been used in a variety of problems. In the application of 
this theory for the temporal stability approach, the frequency was 
considered a constant complex quantity, and the wavenumber k a real 
number. Later it will be shown that, when the imaginary part of the 
frequency is greater tham zero, the disturbance grows exponentially in 
time and is termed a "temporally unstable disturbance." This disturbance 
when amplified to a certain degree gives a definite whistle at a frequency 
which is linearly proportional to the true airspeed. The initial effort 
was devoted to the study of the disturbance caused by sudden enlargement 
in the rear end of the vortex tube as shown in figure 1. 

For the flow described above, the continuity and momentvim equations 
expressed in cylindrical coordinates cire as follows; 

continuity : 


1 _ 3 _ 

r 3r 


(ru) 


r 30 


+ 



( 1 ) 


r-momentum; 


3u . 3u , V 3u ^ 3u v^ 1 3p 

3t 3r r 30 3z r p 3r 


( 2 ) 


0 -momentum: 


3v , 3v . V 3v . 3v uv 1 ^ (3) 

3t " 3r r 30 ” 3z r “ pr 30 


z-momentum; 


3w 3w . V 3w , 3w 1 3p 
3t " 3r r 30 " 3z " ■ p 3z 


( 4 ) 
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tan- 


where u, v and w denote the actual fluid radial velocity V^, 

gential velocity V-., and axial velocity V , respectively. 

y z 

It is the purpose here to derive the equations that control the small 

disturbances of a swirling and steady mean flow. This is accomplished in 

three steps; separation of fluctuations, linearization, and recourse to 

complex functions. 

According to the first step the velocities are separated into a mean 
value and a fluctuation from the mean. It is assumed that in the absence 
of fluctuations the flow is laminar and can be specified as follows: 

u = 0, v = V(r), w = W(r) , p = P(r) (5) 

where the r-component mean velocity is assumed to be zero and the fluid 
beyond the swirling flow generator consists of a forced vortex flow or 
solid body rotation. The mean flow in the 0-direction, the axial mean 
flow, and the pressure terms are assxmed to be functions of radius only. 

If u', v' and w' denote the fluctuations in the r-, 9-, and z- 

directions, respectively, and p' denotes the fluctuation in the pres- 
sure, then the following relationships can be used to express the 
instantaneous velocity canponents and pressure: 

u = u' (r, 9, z, t) 

V = V(r) + v' (r, 9, z, t) (7) 

w = W(r) + w' (r, 0, z, t) (8) 

p = P(r) + p' (r, 9, z, t) (8) 

It should be noted that the mean flow, i.e., V, W, and P, must satify 

the governing equations (1) through (4). Therefore, if equations (6) 
through (9) are introduced into equations (1) through (4) , and the mean 
flow subtracted from these, one obtains; 
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the continuity disturbance equation, 


3u* . . 1 3v' . 3w* - 

9r r r 36 3z ~ 


(10) 


the r-momentvmi disturbance equation, 


3u' 

3t 


+ 


u' 


3u' 

3r 



3u» 

36 


+ 


V' 


iHl) X fw — + w' ^ 
36 / • V” 3z ^ " 3z 


3z 


3z 


+ w' 


3w*\ 
3z / 


2 


Vv' 

r- 


^ _ i 9p' 

r p 3r 


( 11 ) 


the 6-moment\im disturbcince equation. 


3V 

3t 


+ 


u' 


+ w' 


3r 


+ u 


, 3v* 
3r 


3z 


TT- + W' 


3v' 

3z 


— . Zl ^ X — X M 

r 36 r 36 r 36 ” 3z 

Vu' u'v' _ 3p* 

r r ~ ” pr 36 


( 12 ) 


and the z-momentum disturbance equation. 


3w' 

3t 


+ u' 


3r 


+ u' 


3w' V 3w* v' 3w V* 3w' — 3w* 

3r r 36 r 36 r 36 3z 


+ w' 


3w , 3w^ _ _ 1 3p' 

3z 3z ~ p 3z 


(13) 


u * 

In equations (10) through (13), the product terms such as u' • — 
correspond to the effect of one fluctuation on another. Retention of 
these product terms prohibits an easy solution to the equations. However, 
these difficulties are eliminated by assuming that the fluctuations and 
their derivatives have small amplitudes with respect to mean flow quanti- 
ties. This, indeed, is a reasoncible assumption and will be invoked. Thus 
all products and powers (higher than the first) of the fluctuations are 
neglected and only those terms which are linear in them are retained. This 
procedure of linearization leads to the following linearized equations: 



r-momentum : 


Su' V 9u» 

3t r 90 


+ W 


9u' 

9z 


+ w 


. Iw 

90 


W _ 1 ^ 

r " p 9r 


0 -momentum; 


9v' 

9t 


+ u' 


9v V 9v* 

9r r 90 


V* 


90 ” 9z 


+ w' 


9V vu^ 
9z r 


1 9p' 
pr 90 


(14) 


(15) 


z-momentxim; 


9w' 

9t 


+ 


u' 


^ V 9w* 
9r r 90 


r 


w ^ 

90 ” 9z 


w 


9w 1 9p* 

9z ~ p 9z 


(16) 


For simplicity, the axial mean flow w is assumed to be a constant, and 
since the mean flow in the 0-direction is a function of radius only, the 
above continuity and momentum disturbance equations can then be written 
as the continuity disturbance equation [eq. (10)], as well as the 
r-momentum disturbance equation# 


9u» 

3t 


+ 


90 ” 9z 


Vu* ^ _ 1 9p* 
r ” p 9r 


(17) 


the 0-momentum disturbcuice equation. 


9v' 

9t 


+ 


u' 


9v ^ V 9v' ^ - 9v* Vu' 

9r r 90 9z r 


_1_ 9p* 
pr 90 


(18) 


and the z-momentum disturbance equation. 


^ 4. Z iSil 4. M ii£l _ 1 9p* 

9t r 90 9z p 9z 


(19) 
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Finally, the linearity of equations (10) and (17) to (19) enables one 
to obtain periodic solutions in terms of complex functions. Thus the 
system of partial differential equations (10) , (17) , (18) , and (19) can 
be converted into ordinary differential equations. 

The various parameters describing the perturbation have a (z, 0, t)- 
dependence assumed to be given by 

^i(kz + n0 - wt) ^ 20 ) 


where k is the wavenumber of the disturbance in the z-direction, n is 
an integer (which can only be positive, zero or negative) , cuid (i) is a 
constant (which Cctn be complex). In this analysis, k is assumed a real 
number and (O a complex nvunber that will be determined later. Further- 
more, if f(r), g(r), h(r), eind • q(r) were to denote the amplitudes of 

the respective perturbations u', v', w' amd p' whose (z, 6, t) -depend- 

ence is governed by equation (20) , then one could expect to find solutions 
of the type: 


and 


u' 


v' 


w' 


= i f(r) 
= g(r) e 
= htr) e 


i(kz + n0 - ut) 
e 

i (kz + n0 - (jJt) 
i (kz + n0 - Wt) 


p' = q(r) 


i(kz + n0 - wt) 
e 


( 21 ) 

( 22 ) 

(23) 


(24) 


The superscript denotes a preliminary complex solution which ultimately 
leads to a real solution. In other words, one can obtain a set of solutions 
which are the complex conjugates of the preliminary solutions. The pure 
real solutions can then be constructed as follows: 


u' = 1/2 

(O' 

+ u*') 

(25) 

V' = 1/2 

(0- 

+ V*') 

(26) 
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w* = 1/2 (w* + w*’) 


(27) 


and 


p* = 1/2 (p’ + p*’) (28) 

where the superscript * denotes a conjugate of the preliminary solution. 

Substituting equations (21) through (24) into equations (10 ) , (17 ) , 
(18) and (19), and cancelling the common expression, exp [i(kz + n6 - wt)], 
one obtains the following: 

the continuity equation. 


f^- + -- g + kh = 0 
dr r r 


(29) 


the r-momentum equation. 


(w-— V-kw)f-2— g = --^ ^ 
\ r ' r p dr 


(30) 


the 6 -momentum equation. 


(f * 7) ^ 7 • “ * ■ 


pr 


(31) 


and the z-momentum equation. 


^n ^ - (0 + k w) 


(0 + kWjh = -— q 


(32) 


Since the flow at the section of radius R is a forced vortex or 
solid body rotation, the mean flow velocity in 9-direction V can be 
written as 


V = 


r w 

o 


(33) 
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where u is the upstream fluid angular velocity. However, with 
o 


^ V 
dr r 


2 0 ) 

o 


( 34 ) 


and the axial mean flow velocity W considered a constant W , equations 
(29) through (32) become: 

continuity: 


dr 


+ 



- g + k 
r ’ 


h 


0 


(35) 


r-momentum: 


(W 


nW - kW)f-2C0 
o c o 


p dr 


(36) 


0- momentum: 

2W f-( 0 )-n(») - k W ) g = - q 

o o c pr 


(37) 


z-momentum: 


(03 - nO) - kW)h = — q 
o c p 


(38) 


It can be shown from equations (35) to (38) that 


f(r) I 


dh _ 2 n ^o h 

dr”o3-n03 - kW r 

o c 


(39) 


where the parameter K, a function of (k, to) , is defined as 

2 f4 - (w - » W - k vij^y 

r L rt O C J 


K 


2 


(03 - n 03 - k W )' 

o c 


(40) 



Substituting equation (39) into equation (36) yields the following 
expression for g: 



7i + 

1 

0 

3 

C 

1 

3 

k W . 
c h 

h‘ 

_ 

\ * ^2) 

2 0 ) 

r 

n 

r 


L k ' 

0 




It is interesting to observe .that equations (39) , (40) and (41) are exactly 
the same as SuzuJci's equations (197), (198) and (198') (ref. 16). 

Substituting equations (39) and (41) into the continuity equation (35) 
produces the following Bessel's equation of order n in cylindrical 
coordinates 


d^h 

dr^ 


+ 


r r 



0 


The general solution of this equation is of the form 


(42) 


h(r) = ciJ (rK) + caY (rK) (43) 

n n 

where ci and ca are arbitrary constants. It should be noted that in 
the above development no restrictions were necessary and therefore the 
analysis should be simpler and less restrictive. 

Because the value Y tends to infinity as r approaches zero, it 

n 

is necessary to force the constant ca to zero. Equation (43) then 
becomes 


h(r) = c.J (rK) 

In 

The wall boundary condition stipulates no radial velocity, i.e., 
f(r) = 0 at r = R 
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Hence from equations (39) and (44) one obtains 



Evaluation of the derivative of the Bessel function 


(45) 


dJ (rK) 
n 

dr 


= K 



j^(rK) 


n 

rK 



(46) 


and inclusion in equation (45) produces 

r 2n (0 1 

K J , (RK) - ^ J (RK) - — „ T-- J' (RK) = 0 (47) 

^ n-1 Rn (u-noo^- kW^)R n J 

From equation (40) , however, it can readily be shown that 

(0 = In + ^ 1 W + k W (48) 

L 1 + K^ (k"2)J ° 

It should be noted that, if k = 0, equation (47) is automatically 
satisfied, and equation (48) becomes 

(0 = n 0) (49) 

o 

This shows that the angular velocity w is linearly related to the constant 
upstream fluid emgular velocity The angular velocity to can be simply 

related to the soiind frequency as ui = 2 it F, where P is the frequency 
response in kHz. The upstream fluid angular velocity to^ is in reality 
another form of the true airspeed which can be derived easily. Therefore, 
equation (49) indicates that the vortex whistle frequency response is 
linearly proportional to the true airspeed for K = 0. 



If, however, k is not equal to zero, and since ci cannot be zero, 
then the bracketed quantity of equation (47) must be zero, or 


2 0 ) 

o 


0) - n 0) - kW 

o c 


By rearranging equation (50) , (O can be expressed as 


(0 = c (0 + k W (51) 

o c 


) J (RK) = 0 (50) 



where 


C 


n + 


RK 

n 


2 

J . (RK) 
n ~ 1 

J (RK) 
n 


1 


(52) 


Equation (51) , therefore, is the desired general solution showing a 

linear relation between the frequency response of the vortex whistle and 

true airspeed. This, however, is only the situation as long as the values 

of the parameters (O , k, W , a, and K are known. 

o c 

Although equation (51) appears simple, the complex values for angular 
velocity w and parameter K must be recognized. For this reason it is 
difficult to obtain a quantitative solution. Fortunately, a numerical 
procedure (described in Appendix A) can be used to obtain the desired 
solution. The derivations of the expression for the upstream fluid 
angular velocity as well as the constant mean axial velocity 

are included in Appendix B. 
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NUMERICAL RESULTS AND DISCUSSION 


Introduction 

It has been suggested and svibsequently theoretically proven that the 
sensor frequency response is linearly proportional to the airspeed under 
various conditions. The different combinations of vortex tubes and 
swirlers used are tabulated in table 1. The results from computer cal- 
culations will be presented later. The quantitative results which are 
derived in Appendix A prove that the analytic theory works well for all 
combinations of vortex tubes and swirlers tested in the current work. 

It should be noted that the procedure of calculation presented applies 
only to those geometries similar to the ones discussed in this study. 

For other designs the procedure of calculation would have to be modified. 
The quantitative results amd discussions of the theoretical work are 
summarized below. 

Relationship Between Flow Rate and Frequency Response 

Basically, the flow rate versus frequency response has the same 
relationship as airspeed versus frequency response, i.e., the flow rate 
is linearly proportional to the frequency response. For the swirlers 
with the same exit angle but different flutes, the theoretical frequency 
responses are the same, as shown in figures 2 to 9. Under the specified 
condition and configuration, the frequency response changes dramatically 
as the cross-sectional size of theprecessional flow region changes. 

For example, the frequency response for sensor 2 is approximately 3.5 
times greater than that of seMor 1. The frequency response is higher 
with larger exit angle swirlers. It is expected that these sensors 
can be utilized as an instrument for measurement of flow rates in the 
industrial and research areas. 

Effect of Exit Angle on Swirler 

In general, under the saime conditions, the higher the exit angle, 
the higher the frequency response will be. This is true because at 
higher exit angles the tangential component of flow is higher. The 
higher exit angles generate higher angular velocities or revolutions 



which in turn produce higher frequencies. A comparison of the results is 
presented in figures 10 to 30. 

Effect of Wavenumber 

The results for sensor 1 with wavenumbers one, three, five and seven 
are shown in figures 12 to 15, respectively. A comparison of the results 
for each shows no difference, which proves that the frequency response does 
not depend on the wavenumber. This also indicates that assximing a wave- 
number of one when computing the results is justified - 

Effect of Flutes on Swirler 

The effect of flutes imposed on the swirler is shown in figures 31 
to 36. As the number of flutes on the swirler is increased, the intcike 
flow decreases, as does the frequency response. For the same frequency, 
the airspeed sensor designed as a 4-flute swirler needs approximately 12 
percent more flow or speed than that of a 2-flute swirler. Although the 
two-flute swirler would be more economical, the decision as to how many 
flutes the swirler should have should be determined from both design and 
performance criteria. 


Effect of Swirler Core Size D 2 

As observed from figures 12, 21, and 27, the frequency response is 
higher in the small core size of swirler. Specifically, the frequency 
response for the 1.588- and 1.524-cm core size swirlers is much lower 
than for the 0.953-cm swirler. (See table 1 for the dimensions of each 
configuration.) In the experiments for the 0.953-cm core dieuneter 
swirler with 3 flutes whose exit angles are each 45 degrees, a whistle 
was generated. As mentioned in reference 41, some swirlers would not 
generate a distinct whistle as expected. The ideal range for the swirler 
core size D 2 when expressed in terms of the ratio of D 2 /D 1 is between 
0.81 and 0.67. 
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Effect of Cross-Sectional Size of 
the Precessional Flow Region 

The cross-sectional size of the precessional flow region proved to be 
a significant peir^eter in the design of the vortex tube. The frequency 
responses were much higher for the smaller diameter precessional flow 
regions. Figures 37 to 39 clearly show that much higher frequency 
responses occurred in sensors 2, 4, and 6 than were experienced in sensors 
1, 3 and 5, respectively. (Again, refer to table 1 and fig. 1 for dimen- 
sions of sensors 1 through 6.) The 0.635-cm diameter precessional flow 
region produced a frequency response approximately 3 times that of the 
0.953-cm one. It should be noted that the diameter size of the preces- 
sional flow region can be adjusted to obtain a desired frequency response 

range . 


Relationship Between the Mean Axial Velocity and Airspeed 

From equation (BIO), Appendix B, it is concluded that the mean axial 
velocity W is linearly related to the airspeed. Theoretical data for 
different indicators are plotted in figures 40 to 45. The theoretical 
results show the mean axial velocity for the higher exit angle swirlers to 
be less than for the smaller exit angle swirlers. Figure 46 reveals the 
frequency response to be linearly proportional to the mean axial velocity. 


Relationship Between Reynolds Number and Strouhal Number 

Both Reynolds number • Re and Strouhal number St are defined in 
terms of the axial velocity W and diameter of precessional flow region 
D 3 . Except for the singular point at zero airspeed, the predicted 
Strouhal numbers remained constant as the airspeed was increased. (Note 
that, for a given exit angle, the Strouhal number is a constant regardless 
of the number of flutes in the swirler.) The data for six different 
indicators are plotted in figures 47 to 52. when the sensors with the 
same swirler exit angles are compared, one finds the Strouhal numbers are 
the same. Therefore it is reasonable to conclude that the Strouhal n\imber 
is dependent only on a, the swirler exit angle. The approximate values 
of Strouhal number for 45®, 60®, and 75® swirlers are 0.51, 0.88, and 
1.88, respectively. 



APPARATUS AND TESTS 


General Description 

In order to vinderstand the operation and performance of the true 
airspeed sensor, visualization, closed-conduit and wind-tunnel tests were 
undertaken. The purpose of the visualization test was to carefully 
examine the fluid flow phenomenon-precession and observe how the flow 
velocity affects the precessional motion. Another objective was to 
observe and locate a suitable position for picking up the signal. The 
closed-conduit test (flow rate versus frequency response) provided 
information as to how the flow rate could be related to the frequency 
response. These sensors can certainly be compared with commercially 
available flow meters. Finally, the wind-tunnel test or actual flight 
test was conducted to supply the actual data for the relationship between 
the true airspeed eind the frequency response. It is anticipated that 
additional research will be undertaken to improve the design of the 
sensor in the latter two stages. A complete description of the sensor 
setup and test procedures follows. 

Visualization 

In the early stages of investigation, a water model apparatus was 
set up for observing the fluid flow phenomenon-precession. The descriptions 
of vortex tube geometry and experimental arrangement are shown in figures 
53 and 54. In "figure 53, sections I and II, which were made of brass, 
were used for generating the swirling flow. Axial flow was injected at 
section I while eight tangential jet flows were brought into the stream 
at section II. Both the axial and tangential flows were regulated by 
separate valves to produce the desired swirling flow pattern. To provide 
the necessary space for the occurrence of precessional flow, section III 
was made of transparent Plexiglas and constructed as a Venturi-like 
converging-diverging nozzle. Small air bubbles were brought into the main 
stream from the tangential flow inlet. The air bubbles collected in the 
lowest pressure region. The movement of the low-pressure air bubble 
stream revealed the precessional flow for some definite flow rate. It was 
clearly shown that, the higher the flow rate, the higher would be the number 
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of revolutions of the air bubble stream about the axis of the tvtbe. For 
investigation cind observation purposes, four rolls of high-speed movies 
and some slides were taken from different angles. Some of the pictures 
are shown in figtires 55 and 56. The pictures clearly demonstrate that 
the precessional flow whirled around as it left the central part or as it 
entered the enlargement (divergement) section. 

Closed-Conduit Test 

There are many different types of gas flow meters currently in use . 
Some of the meters make use of either the vortex shedding principle or 
the vortex precession principle. Such meters are usually used for low 
flow rates. The indicator designed during this project could also be 
considered as a potential flow meter. It was therefore considered 
desireible to investigate the relationships between the flow rate and the 
frequency response. It was also anticipated that the results from this 
test could be compared with data from wind-txannel tests. 

The initial designs for the vortex tubes and swirlers are shown in 
figures 57 to 59. The vortex tube in figure 58 was used to investigate 
the effect of the. enlargement exit section wall and to note whether a 
gradual or abrupt enlargement made any difference. Swirlers with 4 
flutes and with 45“ , 60“ , or 75“ exit angles were considered. The 
vortex tubes were made of Plexiglas while the swirlers were made of 
aluminum. The test apparatus is shown in figure 60. The VOL-O-FLO flow 
meter (product of National Laboratories Instrument Inc., Model 50) with 
0.283 m® (10 ft^) per minute capacity was used for the measurement of 
flow rate. The flow rate could be read and recorded through a voltmeter 
and converted to the equivalent flow rate unit, m /min. The sound or 
whistle signals were picked up by an electric condenser microphone and 
passed through an cunplifier which was connected to either an oscilloscope 
or a frequency counter. The schematic diagram of the signal amplifier 
is shown in figure 61 , and the electric condenser microphone used in this 
test is shown in figure 62. It is expected that the flow rate is also 
linearly proportional to the frequency response. 



Wind-Tunnel Test 


A small but well-designed wind tunnel was available for testing (see 
fig. 63). The test area is 15.24 cm by 15.24 cm (6 in. by 6 in). A honey- 
comb type screen, 2.54 cm by 2.54 cm grid, was placed in front of the 
tunnel in order to assure viniform flow at the test section. The speed was 
measured by a pitot-static tube which was aligned with the entrcince plane 
of the vortex tube being tested. The frequency responses were also picked 
up by the electric condenser microphone. The arrangement is the same used 
in the closed-conduit test. Different combinations of vortex tubes and 
swirlers were arranged for testing and their dimensions are given in tcdsle 
1 . 
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EXPERIMENTAL RESULTS AND DISCUSSION 


The experimental results, obtained from both the closed-conduit 
tests and the wind-t\innel tests, did reveal a linear relationship between 
the true airspeed/flow rate and the frequency response. All swirlers 
used in these tests had 4 flutes with an exit angle of either 45®, 60®, or 
75®. The actual dimension h on the swirler shown in figure B-1 (Appendix 
B) is shown in figure 64. The comparison between the designed (used in 
the numerical calculations) and actual cross-sectional areas for flow 
passage at section A-A shown in figure B-2 is also shown in figure 65. 

Past investigators had difficulty eliminating noise from the signal. 
During these tests it was discovered that the proper positioning of the 
microphone was critical in obtaining the desired frequency response. This 
inevitably resulted in eliminating noise interference and thus produced 
a stcible and clean signal. This breakthrough made measurements possible. 

It was observed that a definite minimum airspeed value had to be reached 
first for the frequency response to be stable. This was because of the 
unstable precessional flow at low airspeeds as well as the low amplitude 
signals. The minimum airspeed value was obtained by either lifting or 
lowering the microphone. It was therefore necessary to maintain this 
microphone position for all measurements, as the signal would have other- 
wise deteriorated. The arrangement and position of the microphone is shown 
in figure 66. The results obtained from both the closed-conduit test and 
wind-tunnel test are shown and discussed as follows. 

Closed-Conduit Test 

The first step was to test the sensor with compressed air. The 
results for sensors 1 through 4 are shown in figures 67 to 70. The flow 
rate was recorded and calculated according to the calibration sheet 
provided with the flow meter used in this test. When calibrated at 
atmospheric pressure, the formula for flow rate is 

Q = Fi • P 
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where Q is the flow rate of the air (ft^/min) , P is the differential 
pressure (in. of water), and Fi is the constant 5.23. The flow rate 
Q is then converted to the metric system units . 

As observed from the figures, the frequency response for each sensor 
with different swirlers was linearly proportional to the flow rate over a 
definite range. The working ranges for each sensor with a special swirler 
are tabulated in table 2. It was noted that sensors 2 and 4 had wider 
ranges of frequency response than did sensors 1 and 3. The minimum flow 

for generating a stable signal was 0.009 m^/min, and this occurred 
in indicators 2 and 4 with an exit swirler angle of 75°. The minimum 
frequency response picked up by the microphone was 0.70 kHz, and this was 
observed in sensor 3, which had an exit swirler angle of either 45° or 
60°. The mciximum flow rate working linear range was noted to be 0.157 
m /min, and this was observed with sensor 3 having an exit swirler angle 
of 45°. As predicted in the numerical calculations, it was found that 
the higher the exit swirler angle, the higher the frequency response that 
occurred. The best performance among those sensors was for the one with 
an exit swirler eingle of 45°. In general, the sensors with the smaller 
cross-sectional precessional flow areas had quicker responses at the lower 
flow rates. This smaller cross-sectional precessional flow area produced 
a more stable precessional flow and a more pronounced whistle. Below the 
minimum flow rate, the signal was random and the noise level was high. 
Beyond the maximum point, the frequency response declined. It will be 
shown later that the working ranges listed in table 2 will indeed cover 
the maximim desired airspeed of 321.9 km/hr. 

Wind-Tunnel Test 

Unlike the closed-conduit test, when a sensor was tested in the wind 
tunnel , a much higher airspeed was needed to generate a stable frequency 
response. Hence a modified model for the vortex tube was made (shown 
in figure 71) . A convergent nozzle was placed in front of the vortex 
tube. This enabled the airspeed to reach a higher value at the entrance 
to the vortex tube. When compared with the original sensor, this modified 
model was able to generate a higher angular velocity and more stable 
precessional flow. Also, a divergent nozzle was placed in the back of the 
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vortex txabe in order to minimize the outside influence of the back flow 
region. Such ein arrangement, another breakthrough, contributed to making 
measurements possible. A pitot-static tube was placed in front of the 
convergent nozzle to allow measurement of speed. It was anticipated that 
a stable frequency response was not -possible until' the minimum. frequency 
response, as noted in the closed-conduit test with a specified sensor, 
was reached. Based on this test, the minimum airspeed to obtain a -stable 
signal was 70 km/hr (43 . 5 mph) . 

Due to wind-ttinnel limitations described in figure 63, only a few 
measurements were made. However, through correlation of flow rates, 
airspeeds, cind frequency responses, additional data can be extrapolated 
to higher airspeeds. The combinations of test results for true airspeed, 
flow rate, and frequency response are shown in figures 72 to 83 for the 
individual sensors with a specified swirler. The procedure to obtain 
these figures is described in Appendix C. 

It should be mentioned that, for a specified sensor cind swirler, the 
working range for the closed-conduit test is also the working range for 
the wind-tunnel test. The frequency response is lineaurly proportional to 
the true airspeed. When the desired maximum airspeed of 321.9 km/hr (200 mph) 
is reached, the frequency response is still within the working range. 

In figrires 72 to 83, the converted,, dimensional, ntmierical results 
cire plotted for comparison. The converted numerical results were obtained 
by multiplying the theoretical results by the area ratio shown in figure 
65. To avoid any misTinderstanding two points should be mentioned. First, 
the numerical results were calculated assuming inviscid flow. In reality, 
the viscous effects were, however, significant. The deviations between the 
theory and experimental data are clearly evident, especially for sensors 
2 and 4. Secondly, the pitot tiibe position for measuring the airspeed in 
the wind-tunnel test was not at the same location used in the nvonerical 
calculation. The actual airspeed at the entrance of the vortex tube was 
1.84 times higher them that at the point of the converging nozzle. The 
experimental results for sensors 1 and 3 compare favorably with the 
mmerical results when one follows the procedures and measurements 
recommended in this section. It was concluded that the analytical and 
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experimental results did prove that the frequency response was linearly 
proportional to the true airspeed, while quantitative agreement between 
theory and experiment was sufficiently good to allow theoretical prediction 
only for sensors 1 and 3. 
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SUMMARY AND CONCLUSION 


An experimental and ainalytical study of a true airspeed sensor was 
undertaken. The results of the theoretical study proved to be successful, 
in that it was determined that the frequency response was linearly 
proportional to the airspeed. The performance was shown to be heavily 
dependent upon the design of the vortex tube and swirler geometries. The 
smaller the cross-sectional area in the precessional flow region, the more 
sensitive was the frequency response. The exit euigle on the swirler also 
played a very important role in the performance of the sensor. The 
larger the exit swirler angle, the greater was the frequency response. It 
was also observed that as Reynolds niomber increased the Strouhal number 
remained constant. Moreover, it was noticed that the Strouhal ntimbers based 
on were dependent only upon the exit swirler angles. 

The experimental results revealed similar characteristics to those 
obtained through analysis. This true airspeed sensor has three unique 
characteristics: flow phenomenon-precession, vortex whistle, and no 

moving parts. It was also observed that the position of the microphone 
was critical. When the microphone was located close to the inner wall, 
a nondesireible reindom signal was created. On the other hand, when the 
microphone was placed too far from the inner wall, the capability of 
picking up a signal was reduced. Thus, the best microphone position was 
achieved by adjusting the position until a stable signal was obtained. By 
such a procedure it became possible to make the necessary measurements. 

The working rainges of the sensors were beyond the maximtan airspeed 
required. 

Due to the complexity of the problem, the analytical results are 
valid only for inviscid fluid flow. When compared with the .experimental 
result, the deviations noted are belived to be mainly due to the viscous 
effects. Under the cirrangement described under "EXPERIMENTAL RESULTS AND 
DISCUSSION," it was found that the experimental results were in agreement 
with the analytical results for sensors 1 and 3. Although for this 
project the minimum airspeed was 70 km/hr (43.5 mph) , the minimum airspeed 
achievable could be lowered by placing a modified convergent nozzle in 
front of the vortex tube. By such a modification a stable signal could 
be detected even at lower airspeeds. 


In conclusion/ it was found that this true airspeed sensor has a 
linear relationship between frequency response and airspeed and is also 
independent of altitude. The sensor has proven to be relieible, capable 
of high accuracy, and one that could be easily fabricated and maintained. 
Such an insturment has proven to be desiredsle in measuring true airspeeds 
cuid flow rates. 
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APPENDIX A 


NUMERICAL ANALYSIS AND SOLUTIONS 

Since equation (51) contains the precessional angular velocity U), 
a complex ntomber, cind since the arguments of the Bessel functions are 
also complex, a qucuititative solution of equation (51) is difficult to 
obtain. Therefore, cm alternate method will be used to obtain the 
frequency response for a specified airspeed. 

Throughout the following emalysis all but one pcirameter will be 

assumed constant. Let smd respectively, represent the real 

and imaginary parts of the precessional angular velocity such that 

w = oj + ioj. . When svibstituted into equation (20) , one obtains 
r 1 

i (kz + n0 - ut) ^i^ i (kz + n0 - (i) t) . 

e = e • e r - (Al) 

As evident, for values greater than zero, the perturbation amplitude 

increases as increases. The objective of this analysis is to find 

the maximum value of for a specified velocity such that the cor- 

responding 0)^ will be the precessional angular velocity (or frequency) 
response desired. Ecjuation (50) can be written in the following form: 


J , (RK) 
n ^ 1 

(RK) 

n 



(A2) 


In equation (A2) , n will be assigned the value one. The justification 
for this assumption is based on observations in the present work emd 
results of Chcmaud's work (ref. 17). Chanaud tested a vortex tube sub- 
merged in water. He injected air into the vortex tube to observe the 
air bxjbble flow pattern throughout. He noted the formation of a bubble 
line and that the line did not assume a helical form. This can be 
interpreted to imply that the perturbations are constant along a simple 
helix, which would be the case for n = 1 in equation (A2) . Thus, provided 
parameters CO^, W^, k and R are known, equation (A2) becomes 
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J (RK) , 

-2 — |l + 

Jl (RK) RK 


2 0 ) 


(0 - W • - k W 
o c> 


(A3) 


However, from equation (40) with n = 1, K can be written as 


K = 


I ; 

kJ4 0)^ - (w - 0) - k w^) 

\ o o c_ 

Oi - (I) - k W 

o c 


(A4) 


Substituting equation (a 4) into equation (A3 ) , the right-hand side (RHS) 
of equation (A3) becomes 


RHS = 


01-0) - k W 

1 o c 


RK 


2 b) 


V 


1 + 


4 0) - (03 -0) -kw) 

o o c 


2 V 03 - 03 - k W ^ 

\ o c> 


or 


RHS = — 


1 

RK 


03 - 03 - k W 

o c 


V 


4a)^-(a)-w - kw) 

o o c 


(A5) 


The left-hand side (LHS) of equation (A3) is 


LHS = 


J (RK) 
o 

Jl (FK) 


and since 


J (x) 
o 


= 1 - 


2^(11)^ 2®(31)^ 


+ - 


(A6) 


and 


Jl (X) = - 


2^i:.2I 2®.:. 3: 2 ^. 3!. 4! 


+ - 


(A7) 
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the left-hand side of equation (A3) can be expressed approximately as 


1 _ + (RK)** _ (RK)^ 

~ 2^(1!)^ 2‘*(2!)^ 2^(31)^ 

RK _ (RK) ^ ^ (RK) ^ 

^ 2^(11) (2!) 2®(2:)(3I) 


(A8) 


An analysis of the errors due to replacing the exact expression by equation 
(A8) may be found in Appendix D. 

Substituting equation (A4) into eqqation (A8) and simplifying, one 
obtains 


2 

(W - 0) - kw )® - rA-...2 

o c 4 


4 4 

X(0J-£0 - kw)'* + 

o c 64 


X (^ . 0) - kW )^ - 

o c 2304 




0) - kW 

o c 


14 (0^ .(0) - CO - kW 

U o o c J 

[ 2 

4 co^ - (CO - (0 - kW )^’| 

. o o c J 

[4 <4^ - (0. - 0.^ - kW^)“]’| 
nl/2 

) ^ 1 (0) - U) - kW ) ^ 

J o c 


16 


384 



,n3/2 


- (0) - 0) 

- kW 

(w 

o 

c J 




- (w - w 

- kW 1 

(0) 

o 

c J 


o c 


By introducing the following coefficients: 
C - (a + kW 


o o c 

Cl = 4 0)^ - 

o o 

C2 = RK/2 

C3 - R^K^/4 = cl 


(A9) 


(AlO) 
(cont ’d) 
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Cl* = R^kV 16 = c|/2 

C 5 = r'*k'*/64 = C 2/4 

Ce = R®K®/384 = cl/12 

C 7 = R®K®/2304 = c|/36 

Ca = 0) - kW 

o c 


equations (A9) cind (A5) become, respectively. 


LHS = { (U - C ) *’ - C 3 (Cl + 2C £0 - ( 0 ^) (W - C ) 
I o o o 


+ Cs(Ci + 2C W - (03 - C )^ - C 7 (Ci + 2C to - 

o 00 




1/2 


T (C 2(10 - C )MCi + 2C to - ( 0 =^) - C 4 (t 0 - c 

■00 o 

, 3/2 

X (Cl + 2C to - to^) + Ce (w - c ) (Cl + 2C to - to^) 
o 00 


5/2 i 


0 ) + Ce 


RHS 


2C2 


(Cl + 2C 0 ) - 03^) 

O 


1/2 


By equating equations (All) and (A12) and multiplying both sides by 
9 1/2 

(Cl + 2C 0 ) - 0 ) ) / one obtains 

o 


^ (0. 0.) . I 


= { (03 - C ) - C 3 (Cl + 2C 03 - 03 ) (03 - C ) 

o o o 


+ Cs(Ci + 2C 03 - 03^) (03 - C ) 

o o 


- C 7 (Ci + 2C^ 03 - 03^) )t|C2(03 - C^) 


- C4(03 - C (Cl + 2C 03 - 03^) 
o o 


2x2 


+ Ce (to - c ) (Cl + 2C to - to-^) 
o o 


(AlO) 

(concl'd) 


(All) 


(A12) 


(A13) 
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Sin^lifying emd rearreinging according to the power of (U, one obtains 


P(u) = C39 w® + Ci*o + C41 0 )** + C42 + C43 ( 0 ^ 


+ C44 OJ + C45 “ 0 


(A 14 ) 


where the intermediate C constants in the development eure as follows: 

C9 = 9 C| - Cl 

= 3C Cl - 7 

o o 

Cll = 2c; - 3C^Ci 

C 12 = 8C^ - 2Ci 

Cl 3 = 4C Cl - 4C^ + 2Ci 

o o 

Cl 4 = Cl - 4C^Ci 

Cl 5 = 14C^ - Cl 

O 

C 16 = 4C Cl - 16C® 

o o 

Ci7 = sc** - ec^ci 

o o 

Cl a = 4C^Ci - 2C® 
o o 

Ci9 = 9C^ + 4C - 2C1 

o o 

C 20 = 8C Cl - 8C^ - 4C^ 

o 00 

C21 = 4c® - lOC^Ci + cf 
o o 

C 22 = 4C®Ci - 2C cf 

o o 

C 23 = 12C^ - 3Ci 

o 

C 24 * 12 c Cl - 8cl 

o o 

C 25 = 3cf - 12 C^C 1 
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C2 6 " Cz + Cl* + C$ 


C27 = 1 + C3 + Cs + C7 

C 28 = lOC^Cz + C4C9 + C6Ci2 
O 

C29 = 15 C^ + C3C15 + C5C19 + C7C29 
o 

C30 = 10 C®C 2 + C4C10 + C6C13 + 10 C^C 2 C 8 + C4C8C9 + C6C8C12 
o o 

C31 = 20 C^ - C3C16 “ C5C20 * C7C24 

o 

C32 = 5c'*C2 + C4C11 + C6C14 - 10C^C2C8 + C4C8C10 + CeC 8 Ci 3 
o o 

C33 = C3C17 - ISC'* - C5C21 - C7C2S 
o 

C34 = C^C 4 - CiCe 
o 

C341 = c'*C 2(5C8 - c ) 
o o 

C35 = C4C8C11 + C 6 CeCi 4 

C36 = 6c® - C3C18 - C5C22 + 6c C1C7 
o o 

C37 = C^C4 - ClCs 
O 

C38 = - C** + C^CiC3 - cfcs 

O O 

CC = I/C2 
cc 

C39 = ~2^26 ~ C27 

CC 

Cito - “2^26 (C8 - 5C^) + 6C^C2 7 

C 41 = ^(C28 “ SC^C8C26) “ C29 
CC 

C42 = “2^30 + C3I 
CC 

C43 = ■3^32 + C33 

cc 

C44 = ^(C C1C34 + C341 + C35) + C36 

2 O 

C45 = ^(C C1C8C37 - C®C2C8) + C^C 38 + CiC 7 

2 0 O O 
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Equation (51) has therefore been transformed to equation (A14 ) , a sixth 

order polynomial in (0 having real coefficients C 39 through C 45 . As 

stated earlier, for a specified airspeed, one desires the m aximu m value 

of among the six roots of the equation. The corresponding (0^ is 

the frequency response. The solution to equation (A14) will be obtained 0 

through use of a computer porgram available at NASA/LaRC. The subroutine 

of this program, named "RPOLY," is a simple program developed to determine 

the complex roots of an equation by a three-stage variable-shift technique 

(refs. 42, 43). The distinct advantage of this subroutine is that it does 

not require any initial estimates of the roots. 

The derivation of expressions for the upstream fluid cuigular velocity 
0 ) cind the constant mecin cixial flow velocity W is given in Appendix B. 

The computing procedure emd sample programming are given in Appendix E. 
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APPENDIX B 


DERIVATION OF THE UPSTREAM FLUID ANGULAR VELOCITY 

0) AND THE CONSTANT MEAN AXIAL FLOW VELOCITY W 
° c 

The values of the upstream fluid angular velocity to and the dovm- 

o 

stream constant mean axial flow velocity W are strongly dependent upon 

c 

the design of both the vortex tube and swirler flow generator. The 
^s^ivshion will be based on the parameters listed in table 1; however » 
prior to the derivation each parameter will be discussed separately. 

The combination of vortex tube and swirler used in the analysis is 
shown in figure B.-2. The cross-sectional view of A-A is shown in figure 
B-1. The net fluid passage area as shown in figure B-1 is 

\et = I - ni h ~ (Bl) 


where m is the number of flutes in the swirler, and h is the width of 
each flute. It is apparent that the net area at section B-B is 


B 

net 4 


(B2) 


The fluid velocity or airspeed U leaving section A-A at an angle of 
a, as shown below, is separated into axial and tangential velocities U 
and U^ respectively. 



Due to the complexity in geometry, an accurate analysis is very dif- 
ficult to obtain. Thus, an approximation which is close to reality will be 
derived . 



Flute 





1 ** I — 


Scale I L 1 

0.5" 


Figure B-1. Cross-sectional view of section A-A. 



(jj 


Figure B-2. Cross-sectional view of vortex tube and swirler. 



since the fluid passage or gap between the vortex tube wall and swirler 

is so small, it is assumed here that the net area A at section A-A 

net 

is replaced by a circular ring with the same outside radius rj but dif- 
ferent inside radius, named equivalent radius r , which is larger than 

^ t z 

the original radius rz . This equivalent radius r will be determined 

S / 2 

first and then used later as a lower limit in the calculation of conservation 
of angular momentum at section A-A. Hence this circular ring will not 
include the cross-sectional areas of flutes. Thus, one can find this equiv- 
alent radius r by equating the net area A ^ with this new circular 
ring, area as 


A ^ = ir(rf - r^ ) 
net e , 2 


or 



where defined in equation (Bl) . 

Due to the converging geometry between sections, the fluid angular 
velocity at section A-A is different from that at section B-B. The 
conservation of angular momentum can be expressed as 



e,2 



dr 


(B 4 ) 


where w. and to are constant fluid angular velocities at sections A-A 
A o 

and B-B, respectively. 

In the forced vortex flow the tangential velocity is equal to the 

product of the constant angular velocity and the radius. Since the spacing 

between the core of the swirler and inner wall of the vortex tube is 

extremely small, the radius will be considered to be the average of ri 

and r . Thus, the constant angular velocity can be written as 

e,2 A 
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u. 


0 )_ = 


A n + r 


&t2 


1 

2 


or 


2 U sin a 

to = 

A ri + r 

e,2 


(B5) 


Note, however, that = U sin a. Therefore equation (B4) becomes 


r 


2. u sin g ^2 ^ 


ri + r 


e,2 


6/2 


/■ 

•/n 


0) r dr 
o 


(B6) 


The value of the constant angular velocity section B-B can then 

be shown to be 


0 ) = 
o 


r® - 

2 u sin g ^ e ,2 


R" 


ri - r 


e ,2 


(B7) 


Hence, for a given airspeed U and swirler geometry there is one value for 

u . 

o 

The downstream constant mean aucial velocity is obtained by 

considering the continuity of mass for sections A-A and B-B. The conser- 
vation of mass requires 


U A ^ = W B ^ 
a net c net 


(B8) 


where p is constant, U is U cos a, and W is the constant meein 

a c 

axial velocity. Thus, the constant mean axial velocity can be written 

as 


W =s U cos a _ 
c B 


net 


(B9) 


net 
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APPENDIX C 


THE CONSTRUCTION OF FIGURES 72 THROUGH 83 

The procedure for constructing the scales and data line of figures 72 
through 83 is as follows; 

1. Find the slope m of the data line in closed-conduit test 
results . 

2. (A) If the data line goes through the origin, as in (a) , fig;ire 

C-1, calculate the reciprocal of the slope and obtain the 
scale for flow rate per unit frequency response, i.e. 

(m^/min) /kHz. 

(B) If the data line does not go through the origin, as in (b) , 
figure C-1, calculate the reciprocal of the slope and 
obtain the scale for flow rate per unit frequency response. 
Then use the following form to mark the scale for flow 
rate; 

Q = J F + S (Cl) 

m 

3. Mark the scale of the frequency response F on the left-hand side 
of the figure. Mark also the scale of the true airspeed U on 
the bottom of the figure, as shown in figure C-2. 

4. (A) If the data line mentioned in step 2 goes through the origin, 

mark the scale of the flow rate on the right-hand side of 
the figure with zero m /min at bottom line, as in Ca) , 
figure C-2. 

(B) If the data line mentioned in step 2 does not go through 
the origin, mark the bottom line with the value of S 
and the corresponding scales based on equation C-1, as in 
(b) , figure C-2. 

5. Find the slope of the data line taken from the wind-tunnel test, 
then calculate the reciprocal of the slope to obtain the value of 
airspeed per unit frequency response F^. 
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6. Find the intersection of one kHz and the value from step 5. 

7. Draw the line which connects the origin and the intersection 
found in step 6. This line represents the desired relationship 
among true airspeed, flow rate, and frequency response, as 
shown in figure C-3. 
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u 


(b) 


Figure C-2. The scales for flow rate at right-hand side. 
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' APPENDIX D 


ERROR ANALYSIS 

In order to effectively compare the theoretical, numerical, and 
experimental results, an error analysis was undertaken. Throughout 
this research endeavor and, in particular in the development of the 
theoretical aspects of the problem, it became evident that certain 
simplifications would be necessary. As a result of these simplifications 
errors might have possibly been introduced. 

The deviation of the theoretical from the experimental results, 
due to neglecting viscous effects in the development of the theory, has 
already been discussed (see "EXPERIMENTAL RESULTS AND DISCUSSION") . It 
is believed that the viscous effects caused the most significant deviation 
of theoretical results from the experimental data. In the derivation of 
the upstre^ fluid angular velocity 0)^ and the constcint downstream mean 
cixial velocity W^, as presented in Appendix B, the error introduced is 
considered insignificant. The complexity of the flow pattern between 
sections A-A and B-B rendered a quantitative error analysis virtually 
impossible. The derivations, however, were believed to be close enough 
to reality that errors could be considered extremely small compared with 
the overall performance features . 

The most significant errors were found in the derivation of equation 

(A14) , where the Bessel functions and Ji were evaluated on the 

basis of equations (A6) cind (A7) . It is however, very difficult to solve 

the problem, as both arguments of J and Ji are complex numbers. In 

o 

equation (Ad) , J took into consideration the first four terms of 
o 

equation (A6) , whereas Ji took into consideration the first three terms 
of equation (A7) . 

We will now investigate the errors introduced from the elimination of 
the higher terms in the Bessel functions and Ji . The actual values 

of and Ji will be calculated by calling the subroutine COMBES 

provided by NASA/LaRC. This stibroutine COMBES is the only known source 
available for computing the Bessel functions J„^ 2 ) and ^„(z) for 



complex orders and arguments. Note that the convergence of equations (A6) 
and (A7) is premised on the absolute value of the argument of the Bessel 
functions being less than one. The argximent of the Bessel function will be 
derived first and evaluated before proceeding with the derivation. 

It is necessary to first determine the real and imaginary expressions 
for the argument RK. Rewriting equation (40) gives 

- to - n - kw^)' 

K — ^ 

0) - 0) - kw 
o c 


then the expression for RK becomes 


RK = RK 


2 

03 - (n 03 + kW ) + 2(n 03 + kW ) 03 - 03 " 

o o c o c 


(D2) 


03 


- (n 03 + kW ) 


Let Ai = n 03 + kW , A 2 = 4 0 )^ - Ai and A 3 = RK, and stibstitute into 

o c o 

equation (D2) , which then becomes 


RK = A3 


>/a2 + 


2Ai 03 - 0 ) 


(1) - Ai 


(D3) 


Now assume 

03 = 03 + i 03 . (D4) 

r 1 

where 03 and 03 . are the real and imaginary parts of processional angular 
r 1 

velocity, respectively. 

Now substitute equation (D4) into equation (D5) and let 


A 4 = A 3 /[ (03^ - Ai)^ + 03^] 

(D5) 

As = A 2 + 2Ai 03^ - 03^ + 03^ 

(D 6 ) 

As = 2 Ai 03^- 2 03^03^ 

(D7) 
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A7 = 0) • Ai 

r 


(D9) 


One then obtains 

RK = A4 [(A7 - i CO^iCAs + i As)^^^] (D9) 

1/2 

However, it is desirable to convert the expression (As + i As) 
into a simple X + i Y form, where X cind Y are real ncunbers, so that 
the argument RK can be rearranged into a convenient form for later use. 
Noting that Z*^ can be expressed as 

_ ^ln(Z)*^ _ ^aln(Z) _ ^a[ln(Z)+ i(0i + 2m ir)] (DIO) 

where Z=x+iy, a is either a constant or a complex number, |z| = 

2 2 1/2 

(x + y ) , euid 01 is the principle value of angle defined by 

-1 

01 = tan (y/x) . An appropriate selection of 0i allows _ (Ag + iAg) to 
be written as 

1/2 1/2 1/2 

(As + iAs) = e (In (As + a|) + i 0il (Dll) 

where 0i = tan ^ (As/As), or 

0 , 0 , 

(As + iAs) = (As + As) (cos + i sin -^) (D12) 

Again let 


As = (a| + a|) 

(D13) 

0 0 
Ag = A7 COS + co^ sin 

(D14) 

0 0 
Ai 0 = A7 sin cos 

(D15) 

All * A4 As Ag 

(D16) 

Ai 2 = A4 As Ai 0 

(D17) 
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then equation (D9) becomes 


RK = All + i Axa (D18) 

where the constants Ai i and Ai 2 are real and imaginary parts of Idle 
argument. The subroutine COMBES can then be used to evaluate the actual 
values of Bessel functions' and Jx for each argument of RK. 

Now, substituting equation (D18) into equations (A6) and (A7) , one 
obtains 


J (RK) S 1 - 1/4 (All + i Ai 2)^ 
o 


64 


+ i A 12 ) 


1 

2304 


(All + i A12)® 


(D19) 


and 


Ji(RK) = 1/2 (All 


+ i A 12 ) ” (All + A 12 ) 


1 

384 


(All + i A 12 ) 


(D20) 


Expanding the power series, simplifying and separating the real and imaginary 

parts, J (RK) becomes 
o 

J^(RK) = [1 - 1/4 (Ai? - All) + ^ (Ai!( - 6Aif Ail + Ai^) 

(All - 15Ai5t All + ISAif Ai^ - All)] 

1 ^ 3 

l/2Aii Ai2 + (4Aii Ax2 “ 4Aii A 12 ) ’ 

(6A11 Ai2 - 2 OA 11 All + 6A11A1I)] (D21) 


1 

2304 
+ i (- 
1 

2304 


and Ji (RK) becomes 
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Jl (RK) 


(All - 3Aii All) + — (Ai! - lOAif All 


■ [2 

+ 5Aii Ai 2)J + i Ai 2 - (3Aii Ai 2 - All) 
+ 2 ^ (5Aii Ai 2 - lOAii All + Ai|)j 


( 022 ) 


Let 

Ai 3 “All - Ai2 

Ai 4 “All 6A11 Ai2 - Ai2 

AlS = All- ISAii Ai 2 + ISAii Ai2 - A 12 

Ai 6 = “ J All Ai2 

Ai7 = 4Aii Ai2 - 4Aii A ? 2 

Ai 8 = 6 A 11 Ai 2 - 2 OA 11 Ai 2 + 6 A 11 Ai 2 

Ai 9 = 1 “ A 13/4 + Aiif/64 - Ai s/2304 

A 20 “ Ai 6 + Aiy/64 - Aia/2304 


(023) 

(024) 

(025) 

(026) 

(027) 

(028) 

(029) 

(030) 


then- equation (021) becomes 

J (RK) = Ai 9 + i A 20 (D31) 

o 

Note that Ai 9 and A 20 real constant numbers. 

The Bessel ftinction Ji (RK) is obtained in a similar manner to J (RK) 

o 

by allowing 


A21 “ a!i - 3 Aii Ai2 


(032) 
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A22 = All - lOAii Ai2 + 5Aii Ai2 (D33) 
A23 = 3Afi Ai2 - A?2 (D34) 
A 24 “ 5Aii Ai 2 ~ lOAii Ai2 + Ai2 (D35) 
A 25 = y All ” A 22 /I 6 + A 2 2/384 (D36) 
A 26 = Ai 2 “ A 23 /I 6 + A 2 4/384 (D37) 


then the Bessel function Ji (RK) can be written as 


Ji (RK) = A2s + i A2s (38) 

Note that A 2 S and A 26 are real constant numbers. 

Equations (D31) cind (D38) are the approximate solutions for values of 

Bessel functions J eind Ji . Thus far the expressions for the argument 

o 

RK and approximate solutions for and Ji have been esteJslished. 

The absolute values of the approximate Bessel functions euid Ji 

are obtained from equations and (D31) and (D38) , as respectively, as 

E = |j (RK) I = (Ai9 + a|o)^/^ (D39) 

o ' o 


and 


El = |ji(RK)| » (a|s + aIs)^/^ 


(D40) 


The actual value of J as obtained from the computer output is 

o 


J (RK) = BJRE(l) + i BJIM(l) 
o 


from \diich the absolute value is defined as 


A^ + |Jq(RK)I = {(BJRE(l)]^ + [BJIM(1)]^| 


1/2 


(D41) 
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The proced\xre for determining the approximate and actual values of 

J and Ji is as follows: 
o 

1. Determine the values of and from equation (A14) for 

a specified airspeed. 

2. Substitute O) , 0 ). , and other relevant information into 

r X 

equation (D18) to obtain the values of the real emd imaginary 
peurts of the argument, i.e. An and Ai 2 « 

3 . Input the values of Ai i and Ai 2 into subroutine COMBES to 

determine the values of J and Ji and their absolute values. 

o 

4 . Use the value of U) from step 1 to evaluate Ai 9 and A 2 0 

and find the absolute approximate value for as in equation 

(D39) . 

5. Proceed as instructed in step 4 to calculating A 2 s and A 26 
and obtain the absolute approximate value for Ji as in 
equation (D40) . 
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6. 


Compare the solutions from steps 3, 4/ and 5 with the absolute 

actual values of J and Ji as in equations (D41) and (D42) . 

o 

7. Finally, calculate the error for and Ji by using equations 

(D44) and (D45) , respectively- 

The flow diagram is shown in figure D-1, and the computer program 
and a sample output appeeir at the end of the appendix - 

The computer output shows that the arguments of the Bessel functions 
are indeed less than one for each case shown. Therefore, since the 
procedure for deriving and calculating the numerical solutions satisfies 
this requirement, the solution procedure is correct. Furthermore, the 
error analysis reveals insignificant difference between the actual and 
estimated values of the Bessel functions. The errors based on equations 
(D44) and (D45) were 0.11 ^ 10~^ when had the value of 1, and 0.25 x 

10~® when Ji had the value of 0.676 x lo ^ . 
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Figure D-1. Flow diagram for error analysis in the Bessel functions 



COMPUTER PROGRAM FOR DETERMINING VALUES OF J AND Ji 

o 


1 


2 

3 

4 

5 

6 

7 

8 

9 


10 


PROGRAM VmiSTL (INPUT, OUTPUT, TAPE5*INPUT,TAPE6*0UTPUT) 
DIMENSION VEL(ll) ,WC(11) ,RE(11) ,ST(H) ,UU(11) ,WCK(11) , 
DIMENSION U(ll) 

DIMENSION ALPHAD(5) »N1(5) ,PREQ(11) ,XO(ll) 

DIMENSION C0EPF(7) ,R00T1 (6,2) ,Q(11) ,QQ(11) 

DIMENSION BJRE(50) ,BJIM(50) ,YRE(10) ,YIM(10) 

COMPLEX R00T(6) 

NN=11 


CK*1 

MM*6 

WRITE(MM,1) 
FORMAT (IHl) 
ID'EGRE=6 
PI=3. 1415927 
xnu=»i6o.oe-o6 
N»1 


LL=1 

D1=0. 745 / 12 . 

DO 25 LM=1,3 
IP(LM-2) 2,3.4 
D2=0. 625 / 12 . 

GO TO 5 
D2=0.6/12. 

GO TO 5 
D2=0. 375 / 12 . 
CONTINUE 
DO 24 NM=1,2 
IF(NM-l) 6,6,7 
D3=0. 375/12. 

GO TO 8 
D3=0. 25 / 12 . 


CONTINUE 
WRITE (MM, 9) LL 

FORMAT (5X,*FOR THE DEVICE *,Il/) 

R«D3/2. 

H*0.l/12. 

CKW=CK/(2.*PI) 

DO 23 IJ*1,3 

XIJ*IJ+2 

DO 22 11*1,3 

N1(II)*II+1 

ALPKAD(IJ)*XIJ*15. 

WRITE(MM,10) N1(II),ALPHAD(IJ) 

FORMAT (5X, II ,*-* ,F3.,0. ♦DEGREE SWIRLER*//) 
ALPHA*ALPHAD(IJ)*PI/lSO. ^ 

AREAA*PI/4.*(D1**2.-D2**2. )-Nl(n)*H*(Dl-D2)/2. 
ARSAE*FI/4.*D3**2 ^ _ 

RE2= SORT ( ( D1 /2 . ) ♦♦ 2 - ( A REA A /F I ) ) 


(cont'd) 



COMPUTER PROGRAM (Cont'd) 


WRITE(MM,11) ^ .1^ 

11 FORMAT (5X,*VEL*,1 IX, *RE J*,11X,*IM J*,8X,*ERR0R ^*,l4X 
1,*ARGUMENT OF J*/) 

U(1)=0.0 

Q*»»»***4HS>* CHANGE VELOCITY ♦*****♦***♦**♦♦♦*♦♦♦*♦ 

DO 21 J=l,ll 
VEL(J)*U(J)*5280./3600. 

UU(J)=U(J)*1. 60935 
q(j)*vel(j)*gos(alpha)*areaa*6o. 

QQ(J)=UU(J)*COS (ALPHA) *1000. /60.*AREAA*0. 0929 
WC ( J ) * VEL ( J ) *COS ( A LPHA ) *A REAA /A REA B 
WCK(J)*WC(J)*0.304801E-03*3600. 

RE(J)=WC(J)*D3/XNU 

X0(J)=2.*VEL(J)*SIN(ALPHA)*((D1/2,)**3.-RE2**3. )/(R**3 
l.*(Dl/2.+RE2)) 

CO=XO(J)+GK*WC(J) 

C1=4.*X0(J)**2.-C0**2. 

C2=R*GK/2. • 

C3=C2**2. 

C4=G2**3./2. 

G5=C2**4./4. 

C6=G2**5./12. 

C7=G2**6./36, 

C8=X0(J)-GK*WC(J) 

C9=9.*C0**2.-C1 

C10=3>G0*C1-7.*C0**3. 

C11=2.*C0**4.-3.*C0**2?C1 

G12=8.*C0**2.-2.*C1 

Cl 3=4 . *C0*Cl-4 . *C0**3 . +2 . *G1 

C14=C1**2.-4.*C0**2.*G1 

C15=14.*C0**2.-C1 

Cl6=4.*C0*Cl-l6. *00**3. 

C17=9.*C0**4.-6.*C0**2.*C1 
C18=4. *00**3. *Cl-2. *00**5. 

019=9. *00**2. +4. *00-2. *01 
020=8 . *00*01-8 . *00**2 . -4 . *00**3 . 

021=4. *00**3 . -10 .*00**2. *01+01**2. 

022=4. *00**3 . *01-2. *00*01**2 . 

023 = 12 . * 00 ** 2 . - 3. *01 
024=12. *00*01-8. *00**3. 

025=3. *01**2. -12. *C1*G0**2. 

026=02+04+06 

027=1. +03+05+07 

028=10 . *00**2 . *C2+C4*C9+C6*C12 

C29=15.*C0**2. +03*01 5+05*019+07*023 

030=- 10. *00**3. *02+0ii*01 0+06*01 3+ 10. *00**2. *02*08+0^* 

108*09+06*08*012 

031=20. *00**3. -03*Ol6-G5*020-0?*02ii 



COMPUTER PROGRAM (Cont'd) 


CJ2-5 . *C0**4 . *C2+C4*C1 1+C6*C14-10 . *C0**3 . *C2*C8+C4*C8* 
1C10+C6»C8»C13 

033=03*017-15. *C0**4.-C5*C21-C7*C25 
034=00**2. *C4-C1*C6 
0341=00**4. *02* (5. *08-00) 

035=04*08*01 1+06*08*014 

036*6 . *00**5. -03*018-05*022+6 . *00*01**2*07 

037=00**2. *04-Cl*06 

038=-00**4. +00**2. *01*03-01**2, *05 

CC“1 /C2 

039*< 00*026 )/2. -027 
04o=(OC*C26/2.)*(08-5.*00)+6.*00*027 
04l=C0/2. * (028-5, *00*08*026 ) -029 
042=00*030/2. +031 
043=00*032/2. +033 

044=00/2. * (00*01*034+035+0341 )+036 

045=00/2. ♦ (00*01*08*037-00**5. *02*08 )+C0**2. *038+01**3 

1,*C7 • 

00EFF(1)=039 

00EFF(2)=040 

00EFF(3)=041 

00EFF(4)=042 

00EFF(5)=043 

00EFF(6)=044 

00EFF(7)=045 

OALL RPOLY (IDEGRE.OOEFF.ROOT.IERR) 

DO 12 M=l,6 

ROOTl (M. 1 )=REAL(ROOT (M) ) 

R00T1(M,2)=AIMAG(R00T(M) ) 

OONTINUE 

AMAX=R00T1(1,2) 

iiax=i 

DO 13 MN=2,6 

IF(R00T1(MN,2).LE.AMAX) GO TO 13 

AMAX=R00T1(MN,2) 

irAAX=MN 

OONTINUE 

FREQ(J)=ROOTl(IMAX,l)/(2.*PI)/lOOO. 

NOW OALOULATE THE REAL AND IMAGINARY PART OP ARGUMENT 

OP THE BESSEL PUNOTION 

ZR=R00T1(IMAX,1) 

ZI=R00T1(IMAX,2) 

A1=N*X0(J)+0K*WG(J) 

A2=4.*X0(J)**2^A1**2 

A3=R*CK 

A4=A3/((ZR-A1)**2+ZI**2) 

A5=A2+2.*A1*ZR-ZR**2+ZI**2 



COMPUTER PROGRAM (Cont'd) 


IP(J-l) 1^.14,15 

14 BJRE(1)=1. 

BJIM(1)=0. 

BJRE(2)sO. 

BJIM(2)*0. 

A 20*0. 

A25«0, 

A26*0. 

PE0*0. 

PE1*0. 

GO TO 16 

15 CONTINUE 

THETA 1*ATAN2 ( A6 , A 5 ) 

THETA=THETAl/2. 

A7=ZR-A1 

A8*(A5**2+A6**2)**0.25 
A9=A7*COS(THETA )+ZI*SIN (THETA ) 

A 10=A 7*SIN (THETA ) - ZI*C0S ( THETA ) 

A11*A4*A8*A9 

A12=A4*A8*A10 

C ■ FIND THE ACTUAL VALUES OP THE BESSEL FUNCTIONS 
CALL C0MBES(A11,A12,0, ,0. , 1 ,BJRE,BJIM,YRE,YIM) 

C FIND THE ESTIMATED VALUES OF BESSEL FUNCTION. JO 
A13=A11**2-A12**2 
A 14=A 1 1**4-6 . *A 1 2* A 1 2** 2+A 1 2**4 
A 1 5=A 1 1**6- 1 5 . *A 1 1**4*A 1 2** 2+ 1 5 . *A 1 1**2*A 1 2**4-A 12**6 
Al6=-0.5*A11*A12 
A 17=4. *A 11** 3*A 12-4. *A11*A 12**3 
A 1 8=6 . *A 1 1** 5*A 1 2- 20 . *A 1 1** 3*A 1 2** 3+6 . *A 1 2** 5*A 1 1 
A 1 9=1 . -A 1 3/4. +A14/64 . -A 1 5/2304 . 
A20*A16+A17/64.-A18/2304. 

C FIND THE ESTIMATED VALUES OF BESSEL FUNCTION J1 
A21*A11**3-3.*A11*A 12**2 

A22*A11**5-10.*A11**3*A12**2+5.*A11*A12**4 
*A 1 l**2*Ji 12— A 1 2*^'^ 

A24*5.‘*A11**4*A 12-10. *A11**2*A 12** 3+A 12**5 
A 25=0 . 5*A 1 1 -A21/16 . +A 22/384 . 

A 26*0 . 5*A 12-A 23/16 . +A 24/384 . 

C ACTUAL VALUE OF JO 

CA 0*SQRT ( B JRE ( 1 ) **2+B JIM ( 1 ) ** 2 ) 

C ESTIMATED VALUE OF JO 
CEO* SQRT ( A 1 9** 2+A 20** 2 ) 

C ACTUAL VALUE OF J1 

CA 1*SQRT ( B JRE ( 2 ) ♦* 2+B JIM ( 2 ) ** 2 ) 

C ESTIMATED VALUE OF J1 
CE1=SQRT (A 25 ** 2 + A 26** 2 ) 

DCC=ABS(CA0-CE0) 
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COMPUTER PROGRAM (Concl’d) 


C ERROR PERCENTAGE OP JO 
PEO*DCO/CA 0*100 . 

DC1*ABS(CA1-CE1) 

C ERROR PERCENTAGE OP J1 
PE1=DC1/CA1*100. 

16 CONTINUE 

WRITE(MM,17) U(J) ,BJRE(1) ,BJIM(1) ,A11 ,A12,A19.A20 ,PEO 
1,BJRE(2) ,BJIM(2) ,A25,A26,PE1 

17 P0RMAT(5X,P^f.0.2X,*(A) *,E13.7,2X,E13.7»15X,*(*,E13.7 
1,* +I *,E13.7,*)*/11X.*(E) *,E13.7,2X,E13.7,2X,E11.5/ 
211X,*(A) *,E13.7,2X,E13.7/11X,*(E) *,E13. 7,2X,E13. 7,2 
3X.E11.5/) 

IP(J-l) 18.18,19 

18 ST(1)*0.0 
GO TO 20 

19 ST(J)=PREQ(J)*1000.*D3/WC(J) 

20 CONTINUE 
U(7+l)=U(J)+20. 

21 CONTINUE 

22 CONTINUE 

23 CONTINUE 
LIr«LI*l 
WRITE (MM, 1) 

24 CONTINUE 

25 CONTINUE 
STOP 
END 
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(Cont 'd) 
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SAMPLE OUTPUT (Concl'd) 


lao. 

. 200 .. 

1-45, 

vet 

0,' 

20 . 

«0. 

60 r 

DO. 

100 . 


Ui 

(E> 

i*i 

IE) 

_«*) 

IE) 


..10032a6E«01- 

.xiiarafE-oi. 

.inOTa^ErOl. 


.129a9B6EH)2 .i030^£^07. 

.r.Ja««392E-01_. . 

_-.3a66392Er01 .2S387E-oa 



•10032a«E«01 .129a9a6E-02 

.10032a4E*01 .1298986E-02 

jJ X 16l?iy5rai__r»i9ii352E::91. 

.111678*E-0l =.58^%392E-01 


...I .2222344E>01 *I -.1167lo7E*boi 
■ .10503E-07 - 


...25380E-08 


It) jlOO? 2 .a«E*fll t.l 298986 Ez 02 _ 

IE I .1003284£*01 .12989B6E-02 ,1050*E-07 

Ul ■ .111678«Ez01 z.58*^392E-0l 

I E » ~. l X 1 6.7 »W3 8J.EZ.0 8. 


I ' .22 223^^E-0i '»t_j^.ll-671.07EfOOI 


DEGREE _3UI8tER 


ap 'j "iN-J "error « ■“ ARC.0HEH7 OR J ... . 

■ Ui TTnnnnnftP.ni a. I .222236A6-01 »l ~il67ib7E»00) 

IE) .1000000E*01. 0. 0.. 

.Ul 0. 0. 

IE). 0. .Oj 0,. - 

l«l .«nnt>^tP* 0 l‘ fl?qa9a6E-02 1 .2222366Er01. ♦I_z.»n67106E*00) 

IE) * .10032896*0l . .1298987E-02 .10306E*^7 . .. 

Ul .lliaraSErOl . r.98AA39lE-oi._ 

. I E I .111 67836 - 01 -.3 8AA 391E-01 .2 93B3£r08 

Ul .i003286E+bi~_.1298986E^02Z-_T.rrrZ__t. ^2223966-01 .♦! -.ii67i066»00l _ 

IE I .100J2«9et.QI *J 29898 ItzO? ».l.fl5.0.Vfc:fiZ 

Ul .1116785E-0l._r.3866391E.r01 - - 

lEI .1116785E-01 z»58663916-01 ,233876-08 - - 


lAl” .1003286E«01 ^,12989866-02 ..2222366E-0l._»l. r.ll67l06£*00l . 

ISI .1003286E*01. . ,12989876-02... ,103056-07 — ._ 

Ul .lU6785E-01_.r*5866391E-01 

IE) ,11167836-01 -.58663916-01 ,253806-08 .... 

Ul .1003286E601 ,12989866r02_ __L_»i2223.66trbjLj:l_r!>X1^7|066tOOI 

IE I ' '.1003286£*01 ,12989876-02 .. ,105036-07 ...^ 

Ul .1116789E-01 -,5«66391£-01 _ 

I6)....lll67836z01 r,^5B.66391tr01 ,.23383EzC8 

Ul .100J286£»oi .. ,12989866-02 _L .22223666-01 *E r,1167106£*001 

lEI .10017896*0 t ' .12989876-02 ,109056-;07 

Ul" .11167836-Ol_ r.5866391Er01 - - - 

IE I ,11167836r01._ r.5866391tr01 ,25.3e3E-09 
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APPENDIX E 


COMPUTING PROCEDURE AND PROGRAMMING 

Equation (A14) can be easily solved by use of the computer, and the 
step by step procedure is given below. 

1. Input the physical parameters Dj , Dz , D 3 , ni , a and calculate 
equations (Bl) , (B2) , (B3) , (B7) , and (B9) . 

2. Set airspeed equal to 0 km/hr and calculate all the polynomial co- 
efficients listed in Appendix A. 

3. Find the roots of the polynomial by calling the subroutine RPOLY. 

4. Separate the real and imaginary parts of the roots. 

5. Find the maximum value of the imaginary part of root which gives 

the maximum amplitude, and the corresponding real part of root 10 ^. 

6 . Calculate the frequency response and Strouhal number st. 

7. Print out the airspeed and frequency response. Then increase the air- 

speed and repeat the procedure until the airspeed reaches 321.9 km/hr 
(or 200 mph) . 

8 . Repeat the above procedure for different vortex tube configurations 
and swirler designs. 

The flow diagram for the above procedure is shown in figure E-1. All 
the data were plotted by using NASA CDC 6000 series CalComp model 765 
plotter. The NASA/LaRC subroutine INFOPLT was used to generate the plot 
for specified information. The flow diagraim used in obtaining the plot 
is shown in figure E-2. The complete computer program is listed after 
figure E-2. The program only involved sensor 4. For the other indicators, 
only Dz and D$ need to be changed. 


Input the designed parameters 
^ 2 * ^ 3 * ^ 


Calculate 


Set U * 


Calculate Q, W^, Re 


Calculate the numerical values 
of coefficients of equation (66) 


Find the roots of equation (66)^ 
and separate into real and imagi- 
nary parts 


Find the maximum value of Ca)^ and 
the corresponding 


Calculate the frequency from co » 
Calculate Strouhal number St, 


Print out U, F» St 


U + 20 


u - 200: 


Figure E-1. Flow diagram for calculating the frequency at various 
true airspeeds. 









STEP 1 


STEP 2 


STEP 3 



Figure E-2. Flow diagram for getting the plots. 




COMPUTER PROGRAM FOR SOLVING EQUATION (A14) 


PROGRAM WHISTL (INPUT, OUTPUT, TAPE5»INPUT,TAPE6=0UTPUT) 
DIMENSION YEL(11),WC(11),RE(11) ,ST(11) ,UU(11) ,WCK(U) 
DIMENSION ALPHAD(5) ,N1 (5) ,PREQ(11) ,X0(11) ,U(11) 
DIMENSION R00T1(6,2) 

DIMENSION COEPP(7) 

COMPLEX ROOT (6) 

C NUMBER OP POINTS TO BE CALCULATED 
NN*11 

C DEPINE WAVENUMBER K 
CK«1. 

C CHANGE THE VALUES OP K 
DO 
MM=6 

WRITE(MM,1) 

1 PORMAT(lHl) 

C DEGREE OP THE POLYNOMIAL 
IIEGRE=6 

C DEPINE THE PARAMETERS POR THE INDICATOR 
01=0.7^5/12, 

D2=0.6/12. 

D3*0. 25/12. 

R=D3/2. 

H=0,l/12. 

PI=3, 1^15927 
XNU=160.0E-06 

C AN INTEGER REPRESENTS THE PERTURBATIONS ARE CONSTANT 
C ALONG A SIMPLE HELIX 
N=1 

DO 23 IJ*1.3 
XIJ*IJ+2 

C DESIGNED NUMBER OF 7LUTES ON THE SWIRLER 
DO 22 11=1,3 
N1(II)=II+1 
CKW=CK/(2.*PI) 

WRITE(I®1,2) CK,CKW 

2 FORMAT (5X,*THE WAVENUMBER K=*,E10.3,10X,*THE WAVELENG 
1TH=k/(2.PI)=*,E10,3////) 

C DESIGNED EXIT ANGLES ON THE SWIRLER 
ALPHAD(IJ)=XIJ*15, 

A LPHA=ALPHAD ( IJ) *PI/1 8 0 . 

AREAA=PI/4.*(D1**2.-D2**2. )-N1 (II)*H*(D1-D2)/2. 
AREAB=PI/4.*D3**2 
C EQ.(B3) 

RE2=SQRT( (Dl/2, )**2-(AREAA/PI) ) 
WRITE(MM,3)D1,D2,D3,H,N1 (II) ,areaa, areas 

3 FORMAT ( 5X , *D 1=* , FI 0 . 5/5X , *D2=* , FI 0 . 5/5X , *D3=* , FIO . 5// 
15X,»H=*,F10.5/5X,»H0W many flutes in the SWIRLER?*, I 3 
2/5X,*NET AREA AT SECTION A-A=* , FIO . ‘■/SX . *NET AREA AT 
3SECTI0N 3-E=*,?10.5//) 
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COMPUTER PROGRAM (Cont'd) 


WRITE(MM,4) ,N,ALPHAD(IJ) 

k FORMAT (5X,*N=*,I2/5X,*EXIT ANGLE OP SWIRIER** ,P5. 1,2X, 
1*DEGREES*^/) 

U(1)*0.0 

Q********** CHANGE VELOCITY ♦*♦*♦**♦♦*♦♦•♦*♦♦♦**♦♦* 

DO 21 J=l,ll 

G VELOCITY FOOT PER SECOND 
VEL(J)*U(J)*5280./3600. 

C VELOCITY KILOMETER PER HOUR 

uu(j)*u(j)*i . 60935 ^ 

C AXIAL MEAN VELOCITY IN Z-DIRECTI0N(FT/SEC) • 
WC(J)*VEL(J)*GOS(ALPHA)*AREAA/AREAB 
C IN KILOMETER PER HOUR 

WCK(J)=WC(J)*0.304801E-03*3600. 

C REYNOLDS NUMBER BASED ON MEAN AXIAL VELOCITY 
RE(J)=WC(J)*D3ANU 

C UPSTREAM ANGULAR VELOCITY AT SECTION AA AS DEFINED. IN 
C £Q(B7) 

XO ( J) = 2 . *VEL( J) *SIN ( ALPHA ) * ( (DI/ 2 . ) **3 . -RE2**3 . ) /( R** 
13.*(Dl/2.+RE2)) 

WRITE(MM.5)U( J) ,VEL(J) ,UU(J) ,WC( J) ,WGK(J) ,RE( J) ,XO(J) 

5 FORMA T(5X,*U=*,F7.1.2X,*MILES PER HOUR*.5X,*OR*.2X.F7 
1.1,2X,«F00T PER SECOND*, 5X,*OR*,2X.F7.1,2X,*KILOMETER 
2 PER H0UR*/5X,*WC=*,FZ1,2X,*F00T PER SECOND* *5X,* OR*. 
32X,F7.1.2X,*KILOMETER PER HOUR*/5X,*REYNOLDS NUMBER^* 
4,P8.1/5X,*ANGULAR VELOCITY AT ENTRANCE** ,F8. l//) 
CO*XO(J)+CK*WC(J) 

C1*4.*X0(J)**2.-C0**2. 

C2=R*CK/2, 

C3=G2**2, 

C4=C2**3./2. 

C5=G2**4./4. 

C6*C2**5./12. 

C7*C2**6./36. 

C8*X0(J)-CK*WC(J) 

C9*9.*CO**2.-Cl 

C10*3.*C0*C1-7.*C0**3. 

C11=2.*C0**4. - 3 . *C0**2. *C1 

C12*8.*C0**2,-2.*C1 

C13=4 . *C0*Cl-4. *CO**3 . +2 . *G1 

C14=C1**2 . -4. *C0**2 fCl 

C15=14.*C0**2.-C1 

Cl6*4. *C0*C1-16 . *CO**3 . 

C17=9.*CO**4.-6.*CO**2.*Cl 

Cl8*4.*CO**3.*Cl-2.*CO**5. 

C19=9.*C0**2.+4.*C0-2.*G1 
C20*8 . *CO*Cl-8 . *C0**2 . -4 . *CO**3 . 


COMPUTER PROGRAM (Cont'd) 


021=4. *G0**3 . - 10 . *00**2 . *01+01**2 . 

022=4. *00**3. *C1-2.*00*01**2. 

023=12. *00**2. -a*01 
024=12. *00*01-8 . *00**3 . 

025=3. *C1**2. -12. *01*00**2. 

026=02+04+06 
027=1. +03+05+C7 

028=10. *00**2. *02+04*09+06*012 
029=15 . *00**2 . +03*015+05*019+07*023 
030=-10 . *00**3 . *02+04*010+06*013+10 . *00**2. *02*08+C4* 
108*09+06*08*012 

031*20 . *00**3 . -03*016-05*020-07*024 
032=5 . *00**4 . *02+04*011+06*014-10 . *00**3 . *02*08+04*08 
1*010+06*08*013 

033=03*017-15. *00**4. -05*021-07*025 
034=00**2. *04-01*06 
0341=00**4. *02* ( 5 . *08-00 ) 

035=04*08*011+06*08*014 

036=6. *00**5. -03*018-05*022+6. *00*01**2*07 

037=00**2. *04-ei*06 

038=-00**4.+00**2. *01*03-01**2. *05 

GC=1#/C2 

039= *( 00*026 )/2 . -027 
040= ( 00*026/2 .)*( 08-5 . *00 )+6 . *00*027 
04l=C0/2. * ( 028-5 . *00*08*026 ) -029 
042=00*030/2. +031 

043=00*032/2. +033 

044=00/2 . ♦ ( 00*01*03+035+0341 ) +036 
045=00/2 . * ( 00*01*08*037-00** 5 . *02*08 )+00**2 . *038+01** 
13. *07 

0 OOEPFIOIENTS OP THE POLYNOMIAL 
00EPP(1)=039 
00EFP(2)=040 
00EPF(3)=04i 
00EFP(4)=042 
00EFP(5)=043 
00EPP(6)=044 
00EFP(7)=045 

WRITE(MM,6) (OOEPP(I) »I=1,7) 

6 FORMAT(5X,7(E12.5,2X)//) 

0 FIND THE OOMPLEX ROOTS OP THE POLYNOMIAL 
OALL RPOLY ( IDEGRE , OOEFF , ROOT, lERR) 

WRITE(r«M,7) lERR, IDEGRE 

7 FORMAT(5X,*IERR=*,I2,5X,I5///) 

WRITE (MM, 8) 

8 FORMAT (1 IX, *RSAL PART OF ROOT*, 3X, *11^0- PART OF ROOT* 
1,11X,*R0CT NUMESR*/) 



CCMPUTER PROGRAM (Concl'd) 


c 


DO 9 I-1,IDEGRE 
9 WRITE(MM,10) R00T(I),I 
10 FORMAT (5X,2E20. 5,120) 

SEPARATE THE REAL AND IMAGINARY PARTS OF ROOT 
DO 11 M=l,6 


ROOTl (M, 1 )*REAL(ROOT (M) ) 
ROOTKM, 2)=AIMAG (ROOT(M) ) 
11 CONTINUE 


C FIND THE MAXIMUM IMAGINARY PART 0F ROOT AMONG THOSE 
C ROOTS 

AMAX=R00T1(1,2) 

IMAX*1 


DO 12 MN*2,6 

IF(R00T1(MN,2) .LE.AMAX) GO TO 12 
AMAS*R00T1(MN,2) 

IMAX=MN 


12 CONTINUE 

C DESIRED FREQUENCY RESPONSE 

FREQ(J)*R00T1(IMAX,1)/(2.*PI)/1000. 
IF(J-l) 13,13,1^^. 

13 RK=0. 

GO TO 15 

14 0MEGA=SQRT(AMAX**2.+R00T1(IMAX,1)**2. ) 
XK1= OMEGA - 1 . *X0 ( J )- CK*WC ( J ) 
XK2*SQRT(4.*X0(J)**2.-(XK1**2.)) 
XK3=CK*XK2 

XK=XK3/6cK1 

RK*R*XK 


15 CONTINUE 

WRITE (MM, 16 ) RK 

l6^FOmT(5X,*THE ARGUMENT OP THE BESSEL FUNCTION **,E12 

ip(J-l) 17,17,18 

17 ST(1)*0.0 
GO TO 19 

C CALCULATE THE CORRESPONDING STROUHAL NUMBER 

18 ST(J)=PREQ(J)*1000.*D3/WC(J) 

19 CONTINUE 

WRITE(^,20) ROOTKIMAX,!) ,R00T1(IMAX,2),FREQ(J),ST(J) 

20 FORMA T(5X,*THE DESIRED SOLUTIONS ARE*/5X,2E20. 5/5X ,*TH 
IE CORRESPONDING FREQUENCY AT THIS VELOCITY IS*,E15,5,2 
2X,*KHZ*/5X,*THE STROUHAL NUMBER IS*, El 5. 5/////) 

U(J+l)=U(J)+20. 

21 CONTINUE 
WRITE(KM,1) 

22 CONTINUE 


23 CONTINUE 
CK=CK+2. 

24 CONTINUE 
STOP 
END 
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Table 1. Different combinations of vortex tube amd swirler 


Sensor No. 


Size 


1 

2 

Di 

cm 

(in.) 

1.892 

(0.745) 

1.892 

(0.745) 

Dz 

cm 

(in.) 

1.588 

(0.625) 

1.588 

(0.625) 

D3 

cm 

(in.) 

0.953 

(0.375) 

0.635 

(0.250) 


3 

4 

5 

6 

1.892 

(0.745) 

1.892 

(0.745) 

1.892 

(0.745) 

1.892 

(0.745) 

1.524 

(0.6) 

1.524 

(0.6) 

0.953 

(0.375) 

0.593 

(0.375) 

0.953 

(0.375) 

0.635 

(0.250) 

0.953 

(0.375) 

0.635 

(0.250) 



Table 2. The working reinges of the sensors with swirlers having four flutes 
and different exit angles. 


Type of 
Combination 

Flow Rate 
Ranges 

Minimum 

Mciximum 


Q Cm^/min) F 

(KHz) Q (m^/min) F 

(KHz) 

Sensor 1 

45® 

0.0304 

0.72 0.1150 

2.98 


60® 

0.0200 

0.72 0.0680 

3.10 


75® 

0.0200 

0.89 0.0550 

3.05 


Sensor 2 

45® 

0.0152 

1.03 

0.0713 

4.75 


O 

O 

VO 

0.0127 

1.23 

0.0698 

7.50 


0 

in 

0.0090 

0.85 

0.0400 

5.20 


Sensor 3 

45® 

0.0330 

0.70 

0.1570 

3.30 


60® 

0.0200 

0.69 

0.0770 

3.10 


75® 

0.0180 

0.74 

0.0590 

3.05 


Sensor 4 

45® 

0.0206 

1.20 

0.1200 

7.00 


60® 

0.0158 

1.30 

0.0702 

6.65 


75® 

0.0090 

0.84 

0.0450 

5.35 
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Figure 1 . 


Cross“Sectional view of true airspeed indicator • 


FREQUENCY F KHZ 



84 


6 



Figure 3. Flow rate versus frequency response for sensor 1 with 
. swirlers having three flutes and various exit angles. 
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Figure 4 



Flow rate versus frequency response for sensor 1 with 
swirlers having four flutes and various exit angles. 
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Figure 



5. Flow rate versus frequency response for sensor 2 with 
swirlers having four flutes and various exit angles . 
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Figure 6. Flow rate versus frequency response for sensor 3 with 
swirlers having four flutes and various exit angles. 
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Figure 8. Plow rate versus frequency response for sensor 5 with 
swirlers having four flutes and various exit angles. 
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Figxire 9. Flow rate versus frequency response for sensor 6 with 
swirlers having four flutes and various exit angles. 
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Figure 10. True airspeed versus frequency response for sensor 1 
with two-flute swirlers and wavenumber k = 1. 
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Figure 11. True airspeed versus frequency response for sensor 1 with 
three-flute swirlers and wavenumber k = 1. 
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Figure 12. True airspeed versus frequency response for sensor 1 
with four-flute swirlers and wavenumber k = 1. 
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Figiire 13. True airspeed versus frequency response for sensor 1 
with fotir-flute swirlers and wavenumber k = 1. 
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Figiire 14. True airspeed versus frequency response for sensor 1 
with four- flute swirlers and wavenumber k = 5. 




Figure 15. True airspeed versus frequency response for sensor 1 
with four-flute swirlers and wavenumber k = 7. 
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Figure 16. True airspeed versus frequency response for sensor 2 

with swirlers having two flutes and various exit angles. 
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Figure 17. True airspeed versus frequency response for sensor 2 with 
swirlers having three flutes and various exit emgles. 
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Figure 18. True airspeed versus frequency response for sensor 2 with 
swirlers having four flutes and various exit angles. 
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Figure 19. True airspeed versus frequency response for sensor 3 with 
swirlers having two flutes and various exit angles . 
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Figtire 20. True airspeed versus frequency response for sensor 3 with 
swirlers having three flutes and various exit angles. 
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Figure 21. True airspeed versus frequency response for sensor 3 with 
swirlers having four flutes and various exit angles. 
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Figure 22. True airspeed versus frequency response for sensor 4 with 
swirlers having two flutes and various exit angles. 
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Figure 24. True airspeed versus frequency response for sensor 4 
swirlers having four flutes and various exit angles. 
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Figure 25, True airspeed versus frequency response for sensor 5 with 
swirlers having two flutes and various exit angles* 
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Figxire 26. True airspeed versus frequency response for sensor 5 with 
swirlers having three flutes and various exit angles. 
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Figure 27. True airspeed versus frequency response for sensor 5 with 
swirlers having four flutes and various exit angles. 
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Figure 28. True airspeed versus frequency response for sensor 6 with 
swirlers having two flutes and various exit angles. 
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Figure 29. True airspeed versus frequency response for sensor 6 with 
swirlers having three flutes and various exit angles. 
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Figure 30. True airspeed versus frequency response for sensor 6 with 
swirlers having four flutes and various exit angles. 




Figure 31. True airspeed versus frequency response for sensor 1 with 
swirlers having 45® exit angle and various flutes. 
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Figure 32. True airspeed versus frequency response for sensor 1 with 
swirlers having 60® exit angle and various flutes. 
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Figure 33. True airspeed versus frequency response for sensor 1 with 
swirlers having 75* exit angle and various flutes. 
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Figure 34. True airspeed versus frequency response for sensor 2 with 
swirlers having 45“ exit angle and veurious flutes. 



18 



Figure 35. True airspeed versus frequency response for sensor 2 with 
swirlers having 60® exit angle and various flutes. 


FREQUENCY KHZ 



Figure 36. True airspeed versus frequency response for sensor 2 with 
swirlers having 75® exit angle and various flutes. 
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Figure 38. Comparison of size variation in the precessional flow region 
D3 with 4 flutes and exit angle 60® swirler. 


120 



TRUE AIRSPEED U KM/H 


Figure 39. Comparison of size variation in the precessional flow region 
D 3 with 4 flutes and exit angle 75“ swirler. 
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Figxire 40. True airspeed versus meem cuclal velocity for sensor 1 with 
swirlers having four flutes and various exit erngles. 
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Figure 41. True airspeed versus mean axial velocity for sensor 2 with 
swirlers having four flutes and various exit angles. 
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Figure 42. True airspeed versus mean axial velocity for sensor 3 with 
swirlers having four flutes and various exit euigles. 
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Figure 43. True airspeed versus mean axial velocity for sensor 4 with 
swirlers having four flutes and various exit euigles. 
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Figure 44. True airspeed versus mecui axial velocity for sensor 5 with 
swirlers having four flutes and various exit angles. 
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Figure 45. True airspeed versus mean axial velocity for sensor 6 with 
swirlers having four flutes and various exit angles. 
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Figure 46. Relationship between mean axial velocity and frequency 
response for sensor 1 with swirlers having four flutes 
and various exit angles. 
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Figure 48. Reynolds number versus Strouhal number for sensor 2. 
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Figure 50. Reynolds number versus Strouhal nvunber for sensor 4. 
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Figure 52. Reynolds nvunber versus Strouhal number for sensor 6. 
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Figure 54. Hater model setup. 
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Figure 58. Dimensions of the vortex tube with eibrupt enlargement in- the rear. 
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Figure 59. Design and dimensions of the swirler. 
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ELECTRICAL CHARACTERISTICS 

Sensitivity I -65dB ± 4dB (OdB REF IV/ubar at 1 KHz) 
Output Impedance j 600 Ohms 
Omni-directional Response 
Supply Voltagei 2.0 V Min., 10 Y Max, 

Nominal Supplyi 4.5 Y 
Current Draini 0.5 mA Max. 


Figure 62. Description of electric condenser microphone made by 
Radio Shack . 
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Test section 


Figure ^3. Wind tunnel (15.24 cm by 15.24 cm). 





Figure 64. Actual flute width h on swirler used in the experimental test. 








Areas 

of 

Swirlers 

Designed 

Area 

A^ cm^ 

Actual 

Area 

A cm^ 
a 

Ratio. 

of 

• 5^75 cm 

^ 5 ® 

0.7096 

0.6532 

0.9205 

60° 

0.7096 

0.5483 

0.7727 

75® 

0.7096 

0.4838 

0.6818 

D2*1.5240 cm 

45° 

0.8325 

0.7648 

0.9186 

60° 

0.8325 

0.6390 

0.7675 

75® 

0.8325 

0.5616 

0.6745 


Figure 65 . Comparison of designed and actual areas of different 
swirlers with four flutes. 
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Figure 67. Experimental data for 
flutes and various ex: 


swirlers having four 
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Figure 68. Experimental data for sensor 2 with swirlers having fo\ir 
flutes and various exit angles. 
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Figure 70. Experimental data for sensor 4 with swirlers having four 
flutes and various exit angles. 






Figure 71. Modified version of vortex tube for wind-tunnel test. 
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Figure 74. Performance of sensor 1 with 
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c+ 

(0 
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having 4 flutes and exit angle 75®. 
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Figure 76. Performance of sensor 2 with swirler having 4 flutes and exit angle 60®. 
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Figure 82. Performance of sensor 4 with swirler having 4 flutes and exit angle 60®. 
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Figure 83. Performance of sensor 4 with swirler having 4 flutes and exit angle 75“ . 
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